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ABSTRACT 

We study the damping from neutral-ion collisions of both incompressible and compressible magnetohydro¬ 
dynamic (MHD) turbulence in partially ionized medium. We start from the linear analysis of MHD waves 
applying both single-fluid and two-fluid treatments. The damping rates derived from the linear analysis are 
then used in determining the damping scales of MHD turbulence. The physical connection between the damp¬ 
ing scale of MHD turbulence and cutoff boundary of linear MHD waves is investigated. Our analytical results 
are shown to be applicable in a variety of partially ionized interstellar medium (ISM) phases and solar chro¬ 
mosphere. As a significant astrophysical utility, we introduce damping effects to propagation of cosmic rays 
in partially ionized ISM. The important role of turbulence damping in both transit-time damping and gyrores- 
onance is identified. 

Subject headings: magnetohydrodynamics (MHD)-turbulence-ISM: cosmic rays 


1. INTRODUCTION 

Partially ionized plasmas are universal in both space and 
astrophysical environments. They All the most volume 
throughout the Galaxy and build up molecular clouds from 
which stars form. The presence of neutrals affects the 
plasma dynamics indirectly, mediated by t he friction between 
ions and neutrals (see studies by e.g. iPiddingt^ 119561 : 
iKulsrud & Pearcell 19691) . 

On the other hand, astrophysical plasmas are generally 
turbulent and magnetized (s ee e.g. [Armstrong et al.l [19951 : 
iChepurnov & Lazarianl201^ . Both incompressible and com¬ 
pressible MHD turbu lence in fully ionized me dium have 
been e alier studied jGoldreich & Sridhafl I1995L hereafter 
GS95 , iLithwick & Goldreichl 1^0 It ICho & LazarianI |200^ 
l2003h . The GS95 model of turbulence and its extensions 
in compressible fluids successfully reproduce the proper¬ 
ties of MHD turbulence and are supported by numerical 


tests dM aron & GoldreichlUOOlblOio & LazarianlLo02l 120031 : 

iKowal & Lazarianll2O10l). 

The turbulence theory for a partially ionized medium 
should be constructed under the consideration of the damp¬ 
ing by neutrals. There are extensive studies in existing 
literature regarding neutr al-ion collisional d amping, using 
either single-fluid (e.g. iBraginskiil 11965 1 : iBalsaral I199 ct 
[K hodachenko et a l.l2004tlPorteza e t alj2007) or two-fluid de¬ 


script i on (e.g. Pudritzlll990t Balsaralll996t Kumar & Robertd 
20031 IZaqarashvilii et alJ 120111: iMouschovias et al.l 120111 : 
Soler et al.ll2013bllah . The single-fluid approach has simpler 


formalism and hence physical transparency, but is restricted 
on the limited wave frequencies lower than the neutral-ion 
collisional frequency. While the two-fluid approach treats 
ion-electron and neutral fluids separately and is valid over a 
broad range of wave frequencies. 

But all these studies are restricted to linear MHD waves. 
The modern advances of MHD turbulence reveal it is a 
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highly non-linear phenomenon, with the energy injected at 
a large scale cascading successively downwards at the tur¬ 
bulence eddy tu rnover rate (see e.g. iCho & LazarianI 120031 : 
ICho et al]l2003h . The effect of neutral-ion collisions on the 
cas cade of compressible MHD t urbulence was investigated 
by ILithwick & Goldreichl (120011) . but they mainly focused 
on a highly-ionized medium with gas pressure overwhelm¬ 
ing magnetic pressure. The role of neutral viscous damping 
in a mostl y neutral gas in turbu lent motions was thoroughly 
studied in iLazarian et al.l (120041). Combining the two damp¬ 
ing effects, a recent work by IXu et al.l (120141) (hereafter Pa¬ 
per I) carried out a detailed study on damping of the Alfven 
component of MHD turbulence in a variety of turbulence 
regimes and demonstrated its significance in explaining the 
observed linewidth differences between neutrals and ions, and 
the potential value in measuring magnetic held in molecu¬ 
lar clouds. In the current work, we extend the theory de¬ 
veloped in Paper I to include compressible MHD turbulence. 
The weak coupling of Alfven to compressible modes was sug¬ 
gested by GS95. The decomposition of MHD turbulence into 
Alfven, fast, and slow modes has been confirmed and realized 
through a statistical procedure in both Eourier and real s paces 
(ICho & Lazarianll2002L 1200^ iKowal & Lazarianll2010|) . On 

this basis, we can deal with the three components of MHD tur¬ 
bulence separately. We take advantage of the well-established 
linear theory on wave propagation and dissipation in studying 
the damping processes of MHD turbulence. 

The output of this investigation has a wide range of astro- 
physical applications, for instance, cosmic ray (CR) propaga¬ 
tion. The propagation of CRs in ISM is usually described 
in terms of diffusi ve transport of CR s in a turbulent mag¬ 
netic held (see e.g. lSchlickeiseill2002h . The well established 
statistics on the scattering of CRs by magnetic irregulari¬ 
ties can successfully explain the high-degree isotropy of the 
comic radiation. But the conventionally used, simple model 
of CR diffusion assuming a universal diffusion coefficient 
faces big challenges in interpreting more and more observa¬ 
tional findings (e.g . [Hunter et alJll99'^ iGuillian et akllTOOTl : 
lAdriani et al Jl201 U lAckermann et al.l 120 12h . With the de¬ 
velopment and acceptance of the modern concept of MHD 
turbulence (see recent reviews by [Brandenburg & LazarianI 
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1201 4t iBeresnvak & LazarianI l2015h . significant refinements 
and modifications have been progressively made in the picture 
of CR propagation, especially in the aspect of quantitatively 
relating the CR diffusion and the properties of the turbulent 
medium through which they pass. 

The consideration of turbulence damping is necessary for a 
proper description of the interac tion between CR particles and 
MHD turbulence. The study in lYan & LazarianI (120021120041 
|2008|) showed the important influence of turbulence proper¬ 
ties and its damping on CR scattering and propagation in fully 
ionized plasma. The ionization fraction possessed by ISM 
varies over a vast expanse, r anging from fully ionized to al¬ 
most entirely neutral media ([SpitzeilTo^lDrainillonl). In 
the presence of neutrals, the damping process in partially ion¬ 
ized medium is completely different from that in fully ionized 
medium. Accordingly, the propagation of CRs in partially 
ionized medium is likely to present distinctive features and 
deserves a special attention. With this motivation, we proceed 
the earlier research in partially ionized medium and aim to 
achieve a realistic model of CR transport in ISM. 

The organization of the paper is as follows. In Section 2, 
we introduce the concept of neutral-ion and ion-neutral de¬ 
coupling. We then provide a linear analysis on MHD waves 
using both two-fluid and single-fluid approaches in Section 3. 
With the damping rates obtained from linear analysis, damp¬ 
ing of MHD turbulence is presented in Section 4. The numer¬ 
ical tests of our analytical results in diverse partially ionized 
ISM conditions and solar chromosphere are given in Section 
5. Section 6 contains an application on CR propagation in par¬ 
tially ionized medium. Section 7 presents further discussions. 
The summary is given in Section 8 . 

2 . COUPLING BETWEEN NEUTRALS AND IONS 

The coupling of neutrals and ions depends on the wave fre¬ 
quency of interest in comparison with the neutral-ion collision 
frequency, 

Vni = IdPi- (1) 

H ere P i is ion density and 7 ^ is the drag coefficient defined 
in MlflO^ . Vni is used when we consider the drag force 
exerted on the neutrals by the ions. It is related with ion- 
neutral collision frequency Vin by 

^niPn — ^inPi: ( 2 ) 

where pn is neutral density. The condition w ~ i^ni and 
oj ^ Vin correspond to neutral-ion decoupling scale 
and ion-neutral decoupling scale respectively. In a pre¬ 
dominantly neutral medium, is much larger than in- 
At greater scales than neutrals and ions are perfectly 

coupled and oscillate together. In this regime the single-fluid 
approximation is adequate to study the behavior and damp¬ 
ing of MHD waves. Below neutrals start to decouple 

from ions. Perturbations of magnetic field cannot be instantly 
transmitted to neutrals, thus Alfven waves are suppressed in 
neutral fluid. 

The expressions of fcdec are different for different wave 
modes, depending on the form of their respective wave fre¬ 
quencies. In the limit case of Iow-/3, where /3 = 2Cg/V^ 
is defined as the ratio of gas and magnetic pressure, both 
Alfven and fast modes propagate with Alfven speed. Their 
neutral-ion decoupling scales are determined by the condi¬ 
tions VAk cos 9 = Vni and VAk = Vni respectively, where 9 


is the wave propagation angle relative to magnetic field. The 
cos 9 appears for Alfven modes because their waves travel 
along magnetic field, while fast waves propagate isotropically 
and their decoupling scales do not depend on propagation di¬ 
rection. The intrinsic nature of neutrals allows them to sup¬ 
port their own sound waves after they decouple from ions. 
Once satisfying the condition Csnk = Vni, neutrals are able 
to develop their sound waves and become independent from 
ions’ modulation. Thus the of slow modes depends on 
OJ of sound waves. So we can summarize in low-/3 plasma, 
neutral-ion decoupling scale sets the largest scale for the ab¬ 
sence of Alfven waves in neutral fluid for Alfven and fast 
modes, but the presence of sound waves in neutrals in the case 
of slow modes. 

In the meantime, ions are still subject frequent collisions 
with surrounding neutrals above Therefore the re¬ 

maining wave motions in ions can be strongly influenced 
by the collisional friction. On the other hand, the single¬ 
fluid approximation can still be used to describe the MHD 
waves since ions maintain the coupling state with neutrals. 
At shorter scales than ions also become ineffectively 

coupled with neutrals. The scale is determined by oj of 
MHD waves in ions alone. Eor the weakly coupled ions and 
neutrals, the assumption of strong coupling breaks down, so 
the two-fluid description with no restriction imposed on wave 
frequency ranges needs to be adopted. In the following of the 
paper, we use ’’strongly” and ’’weakly” coupled regimes by 
meaning k < /cdec.ni and k > fcdec,in respectively. 

Based on above analysis, Table[T]lists the decoupling scales 
for different wave modes and regimes. 


TABLE 1 

Decoupling scales of MHD waves. 



^dec,ni 

^dec,in 

Alfven 

Uni 

Uin 

COS 9 

VAi cos 6 

Fast 

Uni 

Uin 


Va 

VAi 

Slow 

Uni 

Uin 

Csn 

Csi COS 9 


3. MHD WAVES IN PARTIALLY IONIZED PLASMA 

The study of damping process is based on the linear anal¬ 
ysis of MHD waves in partially ionized plasma. We present 
the wave frequencies obtained from both two-fluid and single¬ 
fluid approaches. And we mainly focus on the solutions in a 
low-/3 limit. 

3.1. Two-fluid approach 

Alfven waves In a neutral-dominated medium, the dissi¬ 
pation of Alfven waves can be enhanced by both netral-ion 
collisional damping and the viscosity effects in neutrals. The 
relative importance of the two damping mechanisms can vary 
in different phases of ISM. When the neutral-ion collisional 
damping is much more efficient and viscous damping can be 
safely neglected, Alfven waves have a well-studied dispersion 
relation fsee e.g. lPiddingtonlll956l : iKulsrud & Pearcd[l969t 

iSoler et al.ll2013bl). 

oj'^ x)^'niW^ — fc^ COS^ ~ iVnik^ COS^ = 0, 

(3) 
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where x = Pn/Pi- Under the assumption of weak damping 
jw/l <C |w_r|, the approximate solutions are 


2 e cos^ 0Vl{{^ + x)vli + e cos^ 0Vl) 
e + + e cos2 0vi 

(4a) 

_ fnixe COS^ 0Vl 

2((1-f cos2 ■ 

(4b) 

It is reduced to 


= V^e cos^ 0, 

(5a) 

^nVje COS^ 0 

W/ = „ 

(5b) 

in strongly coupled regime, i.e. w <C r'm, and 


uj\ = vie cos^ 0, 

(6a) 

^in 

a.r = - — 

(6b) 


in weakly coupled regime, i.e. oj ^ 

The assumption of weak damping holds in both strongly 
and weakly coupled regimes, but fails at intermediate wave 
frequencies. The strongly damped region can be confined by 
the condition jw/l = |u;ri|. FromEq. Q and (|6ll, it yields the 
boundary wavenumbers. 


kt 

kc 


Va^u cos6 »’ 

^in 

2 Vai cos 0 


(7a) 

(7b) 


jw/l converges with at kf. Within [A:+, k~], the Alfven 
waves become nonpropagating and have purely imaginary 
wave frequencies. Although the precise calculation of the cut¬ 
off interval of Alfven waves shou ld involve the discriminant 
of Eq. Q (e.g. ISoler et ^1201 3bl) . numerical results confirm 
the boundaries [fc+, k~] obtained from the equality between 
\ujj\ and |wi^| can very well confine the nonpropagating inter¬ 
val with uiR = 0. Therefore we choose this simple derivation 
and use the resulting kf as the limit wavenumbers of the cut¬ 
off region. 

By comparing with Table [T] we find the cutoff boundaries 
are connected with the decoupling scales by 

II ^ ^dec,ni,||; (8^) 

^c,|| (^^) 


showing the cutoff region is slightly smaller than the 

[^dec.nii ^dec.in] range. 

When neutral viscosity dominates over the neutral-ion col- 
lisional damping, the approximate dam ping rate in strongly 
coupled regime is (iLazarian et al.|[2004l : Paper I) 


\uJl\=T^^ =k‘^iyn, (9) 

where i/n = Csn/(nnUnn) is the kinematic viscosity in neu¬ 
trals. By equaling the above Iw/I and |wi^| from Eq. ( l5al i. the 
corresponding cutoff wavenumber is 


kt = 


Va cos 0 


( 10 ) 


The wave frequencies and k^ in weakly coupled regime are 
the same as those in the case of neutral-ion collisional damp¬ 
ing. 

Fast and slow waves We then turn to the magnetoacous¬ 
tic waves. We consider neutral-ion collisional friction as the 
dominant da mping mechanism. W e adopt t he dispersion rela¬ 
tion g iven bv ISoler et alJ (1201 3a^ (see also lZaaarashvili et alJ 
1201 ill . Taking advantage of the weak damping assumption, 
analytic solutions can be obtained in the limit case of strongly 
coupled fluids. 



(c^ + Vl) ± + Ui)2 


— 4csU| cos^ 0 


2r'm[fc2(c2 + U4) - 2a;|i,] 


(11a) 

(lib) 


Here Cs = 1 / c^^Ci + CsAn is the sound speed in strongly 
coupled ions and neutrals. The a bove solutions a re con sistent 
with those in earlier work, e.g. lEerriere et alJ (I1988I) . Dif¬ 
ferent from Alfven waves, the compressibility introduces the 
important parameter /3. Its definition was given in Section |2] 
In a I 0 W -/3 plasma, the solutions in strongly coupled regime 
have simple expressions. 


= V^k^ 


R — '' A'^ 

^nVie 


LOI = 


‘2iy^ 


for fast modes, and 


wIj = Cg/c^ cos^ 0 , 

^nC^ski 


LOI = - 


2r'r7, 


(12a) 

(12b) 

(13a) 

(13b) 


for slow modes. 

In weakly coupled regime, under the low-^ condition, the 
propagating component of wave frequency becomes 


= Vie 


(14) 


for fast modes, and 

cos^ 0 (15) 


for slow modes. And they have the same damping rate as 
Alfven waves at high wave frequencies (Eq. (l6bb). 

If we compare the expressions of w/ of the three modes 
in strongly coupled regime, we find Alfven and fast modes 
both have the damping rate propo rtional to their quadratic 
wave frequencies (Eq. (1501, (fT20 'l. as they can always in¬ 
duce oscillations of ions orthogonal to magnetic field. The 
additional magnetic restoring force exerted on ions results in 
relative drift between ions and neutrals and damping to the 
wave motions. It indicates faster propagating waves are more 
strongly damped. But for slow modes, we see from Eq. (I13bl l 
that there is no damping to purely parallel propagation with 
0 — 0, since in low-/3 plasma, slow modes propagating along 
magnetic field are analogous to sound waves in the absence of 
magnetic field, and hence undamped. 

We then compare the expressions of loj in both strongly 
and weakly coupled regimes. We find the neutral-ion colli¬ 
sions at low wave frequencies are responsible for coupling the 
two fluids together, and more frequent collisions bring weaker 
damping to MHD waves. Differently, the damping rate in the 
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weak coupling regime ( dSbl ) purely depends on ion-neutral 
collisional frequency. The collisional friction plays the role 
of dissipating waves in ions. 

Analogically, the cutoff boundaries can be determined from 
\^i\ = |w/{|. Fast modes have 


kt = 

k~ = 


^in 


(16a) 

(16b) 


which are the same as Eq. (|7]i for Alfven waves except 
are replaced by kf. Slow modes have 


kt 

kc 


2Vni COS 9 
c-s^n sin 6 

^in 

2csi cos 9 ’ 


(17a) 

(17b) 


The relations between k^ and the decoupling scales of fast 
waves (see Table [T]i are similar to those of Alfven waves in 
parallel direction to magnetic field, 

, 2 

kf^ = ~ A^dec.ni; ( 18 a) 

Sn 

kf^ = “ fcdec.in ■ ( 18 b) 

Slow waves have the same k~^^ = ^fcdec.in.n as Alfven waves, 

but the relation between and fcdec.ni is less straightforward. 
We will discuss this character of slow modes in Section|4] 

In fact, in the context of anisotropic MHD turbulence, the 
wave propagation angle 9 appearing in fc^and decoupling 
scales also depends on scales. In Section!?] we will present 
their exact expressions. 


3.2. Single-fluid approach 

Although the single-fluid approach has a restriction on the 
range of wave frequencies, i.e. w < Vin, it has the advantage 
of a simpler mathemat i cal tre atment. We take the dispersion 
relations from iBalsaral ( 119961) under the single-fluid approxi¬ 
mation. We also follow the approximation of negligible ion 
density applied in their work for this subsection. 

The dispersion relation of Alfven waves is 

-f ikVAn cos 9\ - vl^ cos^ 6 » = 0. (19) 


The same notations are used as defined in their paper. 


A=“,t = l:L,L=14L£5£^, 

k Vni 


( 20 ) 


Notice that Eq. ( fT9] ) can be directly obtained from Eq. (O by 
omitting the first term in Eq. (|3]l, which is unimportant at 
low wave frequencies. 

The exact solutions to Eq. ([T9j) are straightforward. 


A = 


—ik ± 

2 


Vau cos 9. 


( 21 ) 


As is shown by iKulsrud & Feared (Il969h . the solutions be¬ 
come purely imaginary when k > 2, which is equivalent to 
/c|| > 2vni/VAn- From Eq. (jTali, we see 2r'„i/VA„ fc+y in 
mostly neutral plasma. It shows the single-fluid approach also 
can capture the lower cutoff boundary. 


At fc < 2, we express A as {wu -f iu}i)/k. Together with 
Eq. (jail, we get 

<^R = ^ ^ ^Ank'^ cos^ 9, (22a) 

1 ~ 1 /^^ C'C\‘Sp' f) 

oji =--kk cos 9VAn = --• ( 22 b) 

In comparison with the solutions to the two-fluid dispersion 
relation in strongly coupled regime (Eq. (|5]l), with <C 4 
at low wave frequencies, Eq. (l22l) coincides with Eq. (|5]l in 
weakly ionized plasma. 

The dispersion relation of compressible modes is (iBalsaral 

[HU, 

A" + fkVAn sec 9X^ - (cl + VjJA^ 

— ikVAn sec 9c‘^X -I- VjnCn COS^ 9 = 0. 


Recall that we derived the analytical solutions to the two-fluid 
dispersion rela tion under the weak damping approximation 
in Section 13.11 To better capture the wave behavior beyond 
the weak damping limit, here we follow a different procedure 
called Newton’s method to obtain the approximate solutions. 
To simplify the problem, we consider a paradigmatic case of a 
low-/3 plasma. The approximate expressions of wave frequen¬ 
cies in some tractable limit cases are provided in AppendixlBl 
We treat the left-hand side of Eq. (1231) as a function / and 
adopt the approximate roots of / = 0 at the extreme /3 0 

as the starting point, 

,12 —ik sec 0 ± \/4 — gec^ 9 , ^ , 

Aq’ =-^- Vau, (24a) 

Ag’"^ = ±c„ cos 9. (24b) 

Aq’^ and Ag’^ correspond to fast and slow waves respectively. 
Based on Aq, a better approximation is given by 


Ai — Ag — 


/(Ao) 

/'(Ao)’ 


(25) 


where /'(Aq) is the derivative of / at Aq. This one-step itera¬ 
tion can already converge to the actual solutions well. But the 
expressions of Ai can be too complicated to be illustrative. 
Therefore, instead of applying Ai (not shown for simplicity), 
we And Ai can be further reduced at low-^ condition and be¬ 
come 


A 

A 


1,2 

1 

3,4 

1 


—ik sec 9 ± sec^ 9 


^An ) 


, ^ cik tan^ 9 cos 9 

±Cn COS 9 -—- 1. 

2VAn 


(26a) 

(26b) 


We keep A]’^ the s ame as Ag’^ since the modification is neg¬ 
ligible. Erom Eq. (I26al) . we see fast waves become nonprop¬ 
agating when k > 2cos0. That i s, k > 2vni/VAn, in accor¬ 
dance with given by Eq. (I16al i. 

When k < 2 cos 9, we rewrite Eq. (l26l l in terms of and 
uJi, and get 


2 4 — sec^ 0, ,2 ,2 

^- Vi^k^, 


4 

kk sec 9 , 


wj = 


-VAn = -' 


VLk^ 


2'L'ri. 


(27a) 

(27b) 
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for fast waves, and 


■’R = Csnk COS e, 


UJI = -- 


cl kk tan^ 9 cos 9 
2Va„ 


„2 l 2 
''n^± 


(28a) 

(28b) 


for slow waves. In weakly ionized medium, they are in good 
agreement with Eq. (fTSI) and (fOT l derived using two-fluid ap- 
p roach. _ 

IZaqarashvili et al.l (1201 2h argued that the appearance of the 
cutoff wavenumber in Alfven waves is not a physical phe¬ 
nomenon under the single-fluid description. Nevertheless, 
above analysis demonstrates in strongly coupled regime, the 
single-fluid approach provides the same analytical wave fre¬ 
quencies as those obtained using two-fluid approach in neutral 
dominated medium. Furthermore, the different analytical pro¬ 
cedure we applied for single-fluid MHD can inform us with 
the lower cutoff boundary directly from the expression of w/j 
in the cases of Alfven and fast waves. It indicates the single¬ 
fluid approach is adequate to deal with turbulence damping 
when damping occurs in strongly coupled regime. 

In addition, we also need to point out, similar to the sim- 
plifle d expression of the two-fluid ujj at low-/3 con dition (Eq. 
ill3bb l. we also see from the single-fluid ujj (Eq. (I28bb l that 
damping is absent in the case of purely parallel propagation of 
slow modes. This conclusion was earlier reache d through nor¬ 
mal m ode analysis in single-fluid approac h bv iForteza et al 


(2007), whereas the two-fluid approach by 
(1201 Ih and the energy equations given by 


Zagarashvih^ljiL 

Braginskiil (Il96-5h 


both show non-zero damping rate for purely parallel propaga¬ 
tion of slow mode s. Ind eed, according to the general expres¬ 
sion of LUj in Eq. (fTTbb . LUj is not zero at 9 = 0. But it does 
not mean the s impli fied expression of ujj at low-/? me dium 
(Eq. (I13bb and (I28bb ) is invalid. Because first, from Eq. (II Ibb 
we see the term with /cy is subdominant and can be safely ne¬ 
glected at low-/? condition; and second, A:|| becomes more and 
more insignificant compared with k± towards smaller scales 
due to the turbulence anisotropy (see next section). The above 
analytical damping rates will be numerically tested in Section 
0 


and k± when applying the scaling relations, but keep in mind 
that instead of ordinary Fourier components of wave vectors, 
they should be treated as inverse of scales measured parallel 
and perpendicular to the local mean magnetic held. 

The scaling relations were first introduced for trans- 
Alfvenic turbulence in GS95, and later generalized for differ¬ 
ent re gimes of MHD turbulence in cluding sub-Alfvenic turbu¬ 
lence (iLazarian & Vishniaclll999t) . With the scale-dependent 
anisotropy of MHD turbulence taken into account, the above 
linear analysis of wave damping can be used to set the damp¬ 
ing scales of MHD turbulence. We next briefly describe 
the cascading of MHD turbulence and list the corresponding 
damping scales derived in Paper I. 


4.1. Alfven modes 

Alfven modes cascade independen tly from compr essible 
modes and have cascading rate tGS95. lLazarianll2006l) 


fc2/3i-i/3y^^ (29a) 

1/k < Ia, (29b) 


for super-Alfvenic turbulence, i.e., Ma = VlIVa > 1. Here 
Ia is the scale where turbulence anisotropy starts to develop 
and increases towards smaller scales. At scales smaller than 
Ia, the parallel and perpendicular components of the wave 
vector with respect to the local direction of magnetic held are 
related by 

fc|| lA\kj_lAf/\ (30) 

Sub-Alfvenic turbulence has Ma < 1, and a cascading rate 
(ILazarian & Vishni^l 19991) 




{kvif 

LOA 


evi 


kA_VA 

VI/I A. = klL-iVLMl 


Itr < 1/fc < L, (31a) 
l/fc</t„(31b) 


where Itr is the transition scale from weak to strong turbu¬ 
lence. Weak turbulence only evolves in perpendicular direc¬ 
tion with k± increasing during cascade, while strong turbu¬ 
lence has the scaling relation 


4. DAMPING OF MHD TURBULENCE 

MHD turbulence is a highly non-linear phenomenon, as a 
cascade of turbulent energy towards the undamped smallest 
scale. The remarkable feature of the GS95 model of tur¬ 
bulence is the hydrodynamic-like mixing motions of mag¬ 
netic held lines in perpendicular direction have the same 
time scale as that of the wave-like motions along held lines. 
That is, Alfvenic turbulence cascades over one wave period, 
k\\VA ^ k±v±, which is called the critical balance condition. 
Combining the Kolmogorov scaling for the hydrodynamic- 

like motions v±(xkj_ , the critical balance condition nat- 

urally leads to the scale-dependent anisotropy, fey ^ kf . It 
cannot be overemphasized that the GS95 laws are true only 
in local reference system in respect with the local config¬ 
uration of magnetic held lines. The importance of the /o- 
ca l notion was not poin t ed ou t by GS95, but demonstrated 
in iLazman & Vishnia (1 19991) and later numerical l y testi ¬ 
fied in Ic to & Vishniac (120001) : iMaron & Goldreichl (120011) : 
ICho et ^ (120021) . This specification on reference system 
makes the original GS95 picture more self-consistent and en¬ 
tails the dependence on scales of turbulence anisotropy. For 
the sake of tradition and uniformity, we use the notations fcy 


fcy - (32) 

The damping scale of turbulence can be approximately de¬ 
termined by setting |a;/| ^ With both damping effects 
included, the general expression of damping scale is 


^dam,|| — 


-K +Sf)+ 




(33a) 

^dam ^dam,|| \ 1Ak^ am, II j l/Zcd am < Ia- (33b) 

for super-Alfvenic turbulence, and 


Mam, II 


+ ^) + \/(^n + 




-4 


^dam — ^dam,|| .Jl + LM^^k 

dam,II; i/fcd am Itr • 


(34) 


for sub-Alfvenic turbulence. In most situations, only the dom¬ 
inant damping effect needs to be considered. In the case 
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when da mping is do minat ed by neutral-ion collisions, from 
Eq. (fSbl) . ( I29bl ). (l30l l. (I31bb and (l32l) . we obtain 


Regarding the relation between the cutoff boundaries and 
decoupling scales, from Eq. dm . dMt and (IWb . we deduce 


/ 2i/„A 2 1 _3 

fcdam=fy^j (35) 

for super-Alfvenic turbulence when < I a, and 

A^dam = (36) 

for sub-Alfvenic turbulence when < Itr- 
By applying the scaling relations Eq. (l30l l and (l32] | to the 
decoupling scales of Alfven modes in Table[T] we find A:dam in 
above equations are directly related with /cdec.ni by 


A^dam — ( ^ ^ A^dec.r 


(37) 


for both super- and sub-Alfvenic turbulence. And A:dec,in takes 
the form 

fcdecdn = (38) 

VAi 

in super-Alfvenic turbulence, and 

ka.c,n=i^)-^L^MX^ (39) 

VAi 


in sub-Alfvenic turbulence. Here we use the perpendicular 
component k± to represent total k, since we assume the de¬ 
coupling and damping scales are sufficiently small where tur¬ 
bulence anisotropy is prominent, which is a common occur¬ 
rence in astrophysical environments. 

When neutral viscous d ampin g play s a mo re important role, 
a combination of Eq. Q, ([2^, ® , (EB and (EB yields 

^dam — ^ , 1/^dam ^At (40) 

for super-Alfvenic turbulence, and 

fcdam = 1/fcdam < kr. (41) 

for sub-Alfvenic turbulence. The assumption made here is 
the mean free path of neutral-neutral collisions In is relatively 
small compared with I a in super-Alfvenic turbulence and L 
in sub-Alfvenic turbulence, which is reasonable in most ISM 
conditions. 

It is instructive to compare the cutoff wave numbers derived 
from linear MHD waves and the damping scales of MHD 
turbulence. Due to the critical balance, \ujr\ = holds 
for Alfvenic turbulence. Consequently, the cutoff condition 
\ujii I = I w/1 of MHD waves is equivalent to the damping con¬ 
dition = Iw/I of MHD turbulence. Notice that the wave 
propagation direction, cos 6 = k^k, follows the scaling re¬ 
lations given by Eq. (l30l) and (l3^ . and is scale-dependent. 
By taking this into account, the corresponding to Eq. (iTal l 
and ( fTOb are fully consistent with expressed by Eq. (iTSl l. 
(|3^ . (l40l) . and (HTt . In Table|2] we summarize the expressions 
for the cutoff boundaries kf of the three turbulence modes by 
taking the scaling relations of MHD turbulence into account. 
’’NI” and ”NV” represent neutral-ion collisional and neutral 
viscous damping. The asterisks indicate coincides with 
fcdam in corresponding situations. 


2 3 

k^ = (^^ A:(iec,ni) (42a) 

Sn 

A:” = 2 “^fcdec,m- (42b) 

in the case of neutral-ion collisional damping. The difference 
between them is independent of turbulence properties. 


TABLE 2 

Cutoff boundaries of MHD modes 



kt 

kc 

Alfven 

Super 

NI 

Eq. OD* 


NV 

Eq. igO)* 

Sub 

NI 

Eq. l|36)* 


NV 

Eq. ED* 

Fast 

Eq. l|16at 

Eq. l|16bt 

Slow 

Super 



Sub 




4.2. Compressible modes 

iCho & Lazarianl(l2003h found fast modes form in both high- 
and low-/? media and follow a cascade similar to the one 
in acoustic turbulence. The cascading rate o f fast modes is 
(ICho & Lazarianir2003l: lYan & Lazariadl2004l) 

Ma> I, (43a) 

Ma < 1, (43b) 

where Vf is the phase speed of fast modes. It takes the form 


r-i = 


, A: , 1 


Vt = ]jl{ci + Vi) + yici + Vif 

in strongly coupled regime, and 


4CsI42 COS^ 9 

(44) 


Vf = ]J\ (4 + ^ + Vir - cos2 0 

(45) 

in weakly coupled regime. The transition of Vf takes place at 
fcdec.in if the turbulence survives damping. East modes cascade 
radially, so the wave propagation direction can be considered 
fixed during the cascade. The interse ction scale between 
with Vf from Eq. (l44l i and |w/| (Eq. (II IbH corresponds to the 
damping scale 


Aldam — ^ 


L 3 


+ Fj - 21//) 

Vf [?nV:|(c2-v7) + e.c2y/ 
+ V2 - 2VJ) 


Vf 


yVl{cl-Vf)+UclVf 
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for super- and sub-Alfvenic turbulence respectively. In low-/? 
plasma, the above expressions can be written as 

I Ma > l,(47a) 

fcdam = \ 2/3 

[ MA<l.(47b) 

By comparing the damping scale with the neutral-ion decou¬ 
pling scale (Table [T]i, we find 

) (^dec,ni^A) > 1; (48a) 

2 

) ' (A:dec,nii)-3M|, Ma < l.(48b) 



Since the terms (fcdec.ni^yt)”^ and (fcdec.ni?^)”^ are usually 
much smaller than 1, fcdam < ^dec.ni is a common occurrence. 
It means the damping of fast modes takes place before neutral- 
ion decoupling. It is worthwhile to mention that since the 
cutoff scale is always small er th an the neutral-ion decoupling 
scale for fast modes (Eq. (IlSal i). the cutoff does not play a 
role in the damping of fast modes. 

Slow modes cascade passively and have the same as 
Alfven modes (GS95). According to the critical balance, the 
damping condition = jw/l is equivalent to k\\VA = |u;/|, 
while the cutoff condition |w/j| = |tLi/| gives k\\Cs = |u;/|. 
In a low-/3 medium, is always larger than jwi^l on a cer¬ 
tain scale. It means the cutoff of slow waves happens before 
the equilibrium between and \uji\ is reached. Therefore 
the damping scale of slow modes depends on the cutoff scale 
l/fc+ (see Tabled for the expressions). On the other hand, 
slow modes dissipate subsequently after Alfven modes are 
damped out, so their damping scale should not be smaller than 
that of Alfven modes. Accordingly, in combination with Eq. 
(O and (l36l l. we have of slow modes as 


fcdam = min 






3 

2 1 _ 
L^Vl 


when 1/fcdam < I A in super-Alfvenic turbulence, and 


(49) 


kd^m = min 




'Ma, 


‘2h'ri 


L^Vr "m: 


(50) 

when 1/fcdam < hr in sub-Alfvenic turbulence. In fact, the 
ratio between the two terms, of slow modes and fcdam of 
Alfven modes in Eq. (|4^ and (l50l l reveals 


Ks 

^dam,A 


= < 


^ ) (^dec,ni,||^A) 


Cn\fVA 


— )(A:dec,ni,ll^)” 


, Ma> 1, 

(51a) 

, Ma<1, 

(51b) 


where A:dec,ni,|| = v^IVa- The terms (fcdec,m.||^A) ^ and 
(fcdec,ni,||-^^)~^ are much smaller than 1. If the (3 value of the 
medium is not extremely small, < fcdam.A usually stands. 


Thus we can safely take fcdam = k'^ for slow modes in most 
cases. 

Slow modes have the sa me re lation between k~ and fcdec.in 
as Alfven modes (see Eq. (I42bl) '). We then examine the rela¬ 
tion between k'^ and fcdec.ni- The neutral-ion decoupling con¬ 
dition Csnk = Vni Can be written as 


c ^2 ^2 

Sn^sn'^ 


2 


while the cutoff condition is (see Eq. d) 


(52) 


Uclkl 


= Csk\\. 


(53) 


Under the condition of strong turbulence anisotropy, i.e. k ^ 
k± and k± ^ fey, we see the cutoff condition can be reached 
on a larger scale than the neutral-ion decoupling scale. Thus 
slow modes fall i nto a diffe rent s ituation from fast and Alfven 
modes (see Eq. (II Sal) and (142ab l and can have < fcdec.ni- 
But in terms of damping, due to fcdam = k'^, slow modes are 
similar to fast modes and also damped out when neutrals and 
ions are strongly coupled together. 

Comparing this section with Section[2 we can see the main 
difference between the damping processes of linear MHD 
waves and MHD turbulence is the consideration of scale- 
dependent anisotropy, which depends on turbulence proper¬ 
ties, e.g. super- or sub-Alfvenic turbulence, and turbulence 
regimes, e.g. weak or strong turbulence. For our study on 
turbulence damping of Alfven and slow modes, a fixed wave 
propagation direction, e.g. purely parallel or perpendicular 
propagation, has no physical meaning. Furthermore, as re¬ 
gards scaling estimations, e.g. evaluation of the damping 
scales, turbulence anisotropy is the essential physical ingre¬ 
dient. Ignorance of it can lead to severely wrong conclusions 
in related astrophysical applications. 


5. NUMERICAL TESTS IN DIFFERENT PHASES OF 
PARTIALLY IONIZED ISM AND SOLAR 
CHROMOSPHERE 

The generality of the two-fluid approximation allows us to 
apply our results to a wide variety of astrophysical situations. 
We next apply the analysis on turbulence damping in different 
phases of partially ionized ISM. To avoid confusion, we stress 
that by numerical tests, we mean numerically solving the full 
dispersion relations with given physical parameters, but not 
performing any MHD simulations. 


TABLE 3 

Parameters used for different phases of partially 

IONIZED ISM AND SOLAR CHROMOSPHERE. 



WNM 

CNM 

MC 

DC 

SC 

niiicm 

0.4 

30 

300 

10* 

4.2 X 10*3 

rielna 

0.1 

10-3 

10-* 

10-® 

1.78 X 10-3 

T[K] 

6000 

100 

20 

10 

6220 

B[^l G] 

8.66 

8.66 

8.66 

86.6 

6.96 X 10’* 

P 

0.22 

0.23 

0.20 

0.03 

0.03 

Ma 

0.4 

2.9 

9.2 

5.3 

0.4 


Table [3 lists the typical physical parameters for warm neu¬ 
tral medium (WNM), cold neutral medium (CNM), molecular 
cloud (MC) and dense core in a molecular cloud (DC). The 
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param eters nn, ne/nn, and T are taken from iLazarian et akl 
( 120041) . We assume temperatures are equal for all species. 
Magnetic field strengths are v^x the half m aximum values 
of the Zeeman measurements from figure 1 in iCrutcher et al.l 
(120101) as the total strength of the three-dimensional mag¬ 
netic field. In addition, we ado pt the drag coeffici ent as 
= 3.5 X 10^^cm^g“^s“^ from iDraine et al.l (ll983lfT. and 
the driving condition of turbulence 

L = 30pc, 14 = 10kms-\ (54) 

The masses of ions and neutrals are rrii = = mn for 

WNM and CNM, and r rii = 29mH, = 2.3mH for the 

other environments (IShulll99^ . 

Not only ISM, observations confirm MH D waves are also 
ubiquitous in the solar atmosphere (e.g. iDe Pontieu et al.l 
1200^ l20l^ . The damping process of MHD waves in the par¬ 
tially ionized solar atmosphere, i.e. photosphere and chromo¬ 
sphere, has been widely suggested as a heating mechanism of 
the external solar atmosphere (e.g. IOsterbrockll961l!Hollwegl 
Il986t iNarain & Ulmschneidenll996t lGoodmanll2000l 1200 ll) T 

Thus we also extend our analysis to the solar chromosphere 
(SC) environment. In reality, the physical parameters vary 
with altitude from the photosphere level. For the sake of sim¬ 
plicity, we neglect the vertical variation, but consider a toy 
model in a homogeneous medium with all parameters taken at 
a median height. This simplified analysis serves as a paradig¬ 
matic case of damping in chromosphere-like environment. A 
more sophisticated model may be required when making com¬ 
parisons with observational data. 

We adopt model C for quiet sun in IVernazza et al.l (Il981h 
at a height 1280 km and only consider hydrogen. Magnetic 
field strength is obtained by assuming B = ilph(p/pph)°'^ 
(iLeake & Arbed 1200^ . With Bph = 1.5 x 10^ G and 
Pph = 2.73 X 10“^ g cm“^ at the photospheric level 
used, we get B = 69.6 G. For the neutral-ion collisional 
cross-sec tion, we use the value (ft! 1 x 10“^"^ cm^) pro¬ 
posed by IVranies & Krstid (l2013h . Their approach with en¬ 
ergy dependence and quantum effects contained yields a col¬ 
lisional cross-section two orders of mag nitude larger tha n 
that obtained from the hard sphere model (iBraginskiill 1 9651) . 
Then drag coefficient can be calcul ated with the expression 
(IBraginskiill 1965trSoler et al.ll2013ah . 


1 IIQUbT 

2mH V Trm/r ' 


(55) 


Besides, we use the height of chromosphere as the injection 
scale L = 2500 km, and assume turbulent velocity at L as 
14 = 30 km s“^ based on the velo city amplitude measure¬ 
ment o f the chromospheric waves by lOkamoto & De PontieuI 
(1201 ll) . The parameters are also listed in Table[3] 

According to the values of /? listed in Table [3] in the fol¬ 
lowing comparisons with numerical results, we will use the 
simplified analytical results introduced in previous sections 
for dealing with low-/3 media. 


5.1. Alfven modes 

We assume the energy is equally distributed into Alfven, 
fast, and slow modes at the scale of turbulence driving. 

^ This constant 7 ^; applies in conditions like molecular clouds where the 
relative velocities between neutrals and ions are relatively low and the ef¬ 
fective cross-sect ion of neu tral-ion collisions is much larger than geometric 
cross section (see lShril992ll . 


Fig. 1 1(a) illustrates the damping of the sub-Alfv enic t urbu- 
lence in WNM. The cascading rate is given by Eq. (EH). No¬ 
tice the scales shown are within the strong turbulence regime. 
We show the damping rates obtained by numerically solving 
the two-fluid dispersion relations with only neutral-ion col¬ 
lisional damping (Eq. ([Ji) and both neutral-ion collisional 
and viscous damping (see Paper I). In this case, neutral-ion 
collisions and neutral viscosity have comparable damping ef¬ 
ficiencies in the strongly coupled regime, so the general ex¬ 
pression of fcdam (Eq . (l34l) ) applies. The analytical damping 
rate is from Eq. (l4bl) . which coincides with the numerical re¬ 
sult at both small and large wavenumbers, but underestimates 
the numerical result at intermediate wavenumbers when jw/l 
and jwfll are of the sa me order. The single-fluid damping rate 
is given by Eq. (I22bl l. which traces the numerical one well at 
small wavenumbers. The numerical solutions do not exhibit a 
nonpr opaga ting interv al with ujr = 0. 

Eig. 1 1(b) and |l(c)| display the damping of super-Alfvenic 
turbulence in CNM and DC. The solutions with both damp¬ 
ing effects are in full coincidence with those considering 
only neutral-ion collisional damping due to negligible viscous 
damping. Accordingly, fcdam for neutral-ion collisional damp¬ 
ing (Eq. (l35]) l applies, which coincides with the lower cutoff 
boundary By comparing with numerical results, we are 
convinced that the single-fluid approach can provide a good 
approximation of the actual damping rate down to the damp¬ 
ing scale (or cutoff boundary) of Alfven modes. The pure 
imaginary solutions are omitted from the numerical results, 
corresp ondin g to the cutoff regions. 

Eig. |l(d)| shows the damping of Alfven modes in sub- 
Alfvenic SC. The filled circles represent the analytical damp¬ 
ing rate including both damping effects (see Paper I), 


, I ^ [t^ ^ -f (1 -f x)Vni) + COS^ 0VA^] X^m 

^ 2[/c2 COS^ evl^ + XTv^Vni + + (1 + X)i^m)^] 

(56) 

Obviously neutral viscosity dominates over neutral-ions col¬ 
lisions in damping Alfven modes. It meets the prediction in 
Paper I that neutral viscosity can play a more important role 
of damping in sub-Alfvenic turbulence than in super-Alfvenic 
turbulence. The exact criteria for determining the relative im¬ 
portance between neutral viscosity and neutral-ion collisions 
in damping super- and sub-Alfvenic turbulence can be found 
in Paper I. As in the case of WNM, the general expression by 
Eq. (l34l l including both damping effects can give a good ap¬ 
proximation of fcdam- But here we present /cdam with a simpler 
form (Eq. (|4TJ) derived from viscous damping alone as an 
approximate /cdam- We see in SC condition, owing to the ex¬ 
istence of neutral viscosity, the damping scale is considerably 
larger than that predicted by neutral-ion collisional damping, 
and is also larger than the neutral-ion decoupling scale. It 
shows not only neutral-ion collisions, neutral viscosity should 
also be considered as an important mechanism in damping 
Alfven modes. 

Besides damping rate, the propagating component of wave 
frequency is also present in Eig. |2]as an example in the 
case of Alfven modes in MC. The same symbols are used as 
in Eig. [U except that the dashed line and open squares are the 
numerical and analytical (Eq. (fSal) and d^ l respectively. 
Notice that as pointed out earlier in Section |4l due to the crit¬ 
ical balance, jwfll in strongly coupled regime coincides with 
the cascading rate. We And Iw/jI starts to decay at A:dec,ni, and 
is cutoff at k^. Then arises again at k~, but is only fully 
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resumed after reaching fcdec.in- Thus the region [A:dec,ni, fcdec.in] 
contains the cutoff zone [fc+, k~]. Different from wave damp¬ 
ing, which has damping rate insensitive to decoupling scales, 
the propagating behavior of waves depends on properties of 
fluid coupling, and can only be studied outside [fcdec.ni, ^dec.in] 
region. 

5.2. Fast modes 

Fig. |3(a)P(d)] display the comparison between the analyt¬ 
ical and numerical damping rates of fast modes in different 
ISM conditions and SC. The cascading rate is from Eq. (Il3]l, 
with the wave propagation angle fixed at 45°. The ana lytica l 
two-fluid (Eq. (112bl i and d&bl i) and single-fluid (Eq. (127bl i') 
damping rates are in excellent agreement with the numerical 
result on scales larger than the cutoff boundary l/fc+. 

Both |wi^| and |w/| are shown in Eig. |4]for fast modes in 
MC. The numerical |u;rj| is obtained by numerically solvin g 
the dispersion relation (equati on ( 51) in ISoler et al.ll201^ . 
Analytical \ojr\ is from Eq. (112ab and (fTO . Similar to the 
case of Alfven modes, the cutoff interval is contained within 
[fcdec.ni, fcdec.in], and Only outside [fcdec.ni, fcdec.hj, fast waves 
have phase speed equal to the Alfven speed. 

In all the conditions present, fast modes are damped out 
when neutrals and ions are strongly coupled, i.e. fcdam < 
fcdec.ni- In particular, fast modes in WNM are the most severely 
damped among the ISM phases. WNM has much lower to¬ 
tal density than other ISM conditions, which leads to a faster 
wave phase speed and lower neutral-ion collisional frequency. 
Accordingly, fast modes in WNM have a lower cascading 
rate (Eq. but a higher damping rate (see discussions 

in Section 13.1b . It means the intersection between the two 
rates can take place at a large scale. In SC, high wave phase 
speed comes from strong magnetic field strength and leads to 
a low cascading rate. Although large ion density yields a high 
neutr al-ion collisional frequency and low damping rate (Eq. 
(I12bb ). it can still exceed the substantially slow cascade on 
a large scale. Therefore fast modes in SC are also severely 
damped. The vertical dashed lines indicate the damping scale 
given by Eq. gTli. It shows the approximate fcdam at low-/? 
limit is consistent with numerical results. In fact, according 
to Eq. gi, fcdam has a dependence on the wave propaga- 
tion d irection. Eig. [Sjpresents the damping scale given by Eq. 
(I46bb as a function of 9 in WNM as an example. The damp¬ 
ing scale in other cases show similar results. We find despite 
the existence of anisotropy, the dependence of fcdam on 6 is so 
weak that a constant damping scale can be applied. It leads to 
a conclusion that the energy distribution of fast modes keeps 
isotropic over all existent scales. 

5.3. Slow modes 

Eig. |6(a)||6(d)1 show the results of slow modes. The ana¬ 
lytical damping rate for two-fluid description is given by Eq. 
(I13bb for k < and Eq. (l6hb f or k > k~, and the single¬ 
fluid dam ping rate is from Eq. (128bb . The analytical |wfl| 
from Eq. (113ab for fc < and Eq. (flSl l for k > k~ is also 
shown in Eig. |7]in the case of MC. Erom Eig. |7] we find sim¬ 
ilar to other wave modes at high wave frequencies, although 
lor reemerges at the cutoff boundary k~, the phase speed can 
only reach Csi at the slightly smaller ion-neutral decoupling 
scale. 

Slow modes exhibit a distinctive feature on wavenumbers 
above fcdec.ni- Besides the slow modes in ions, a new sort of 
slow modes emerges in neutrals, which is the ’’neutral slow 


mode” reported in IZaaarashvili et al.l (1201 ll) . It is also the 
’’neutral acoustic mode” coupl e d with the magnetoacoustic 
waves discussed in ISoler et ^ (l2013al) . It is generated be¬ 
low the neutral-ion decoupling scale when the slow modes on 
larger scales impose compression and produce perturbations 
in neutral fluid. In the rest of the paper, we will term the two 
branches of slow modes below the scale 1 /fcdec.ni as ’’neutral” 
and ”ion” slow modes respectively, and term the slow modes 
in strongly coupled regime and ”ion” slow modes together as 
’’usual” slow modes. We next put our focus on the neutral and 
ion slow modes within the wavenumber range [fcdec.ni, fcdec.in]- 
If we extend the wave frequencies of slow modes in 
strongly coupled regime (Eq. (fT3b ) to the scale fcdec.ni, we find 
they reach the values 

Wfl(fcdec.ni) « COS^ 6, (57a) 

w/(fcdec.ni) ~ sin^ 9, (57b) 

where we use the approximations ^ 1 and Cs„ ^ Cg- Start¬ 
ing from fcdec.ni, the slow modes arising in neutrals have wave 
frequencies as 

Wfl = cLfc^ (58a) 

u>i = —sin^ 9. (58b) 

The propagating component ujr has the expression corre¬ 
sponding to pure acoustic waves in neutrals. The i magi - 
nary part uji has the same value as a;/(fcdec.ni) in Eq. (I57bb . 
By setting \u!r\ = jw/j, we can get its cutoff wavenumber 
Vnij i^Csn), which is equal to fcdec.ni/2. Below fcdec.in, ions 
start to separate from neutrals, w/ of the neutral slow modes 
becomes 


The above analyticallo;/1 (Eq. (I58bl l and (fS^ l are shown by 
filled circles in Eig. |6|and|2] Due to the strong anisotropy at 
relati vely small scales, sin 9^1 and the difference between 
Eq. (I58bl) and ( |59] |) is indistinguishable on the plots. The 
filled squar es in Eig. |7] specifically show the analytical jwijj 
(Eq. (I58al i). Although the slow modes in neutrals have the 
properties of pure acoustic waves, it is induced by the slow 
modes in strongly coupled regime and is one of the solutions 
to the dispersion relation of the magnetoacoustic waves in par¬ 
tially ionized plasma. We treat it as an additional branch of 
slow modes in our analysis. 

The relative damping efficiency of the ion and neutral slow 
modes at high wave frequencies depends on the ionization 
fraction of the plasma. The ratio of their damping rates at 
fc > fcdec.in IS x(= Pn/Pi)- The higher degree of ionization 
in WNM and SC results in a smaller ratio between the two 
damping rates compared with other phases. 

Another difference between WNM and other environments 
is th e abs ence of cutoff region in WNM. Also, we see from 
Eig. |6(a)| that the cascading rate of slow modes is above the 
damping rate of ion slow modes over all the scales. Thus 
the usual slow modes in the case of WNM survive neutral- 
ion collisional damping. Eor the rest conditions, the damping 
scales all coincide with the cutoff boundary (see Table |2] 
for the expression). 

We next turn to the ion slow modes within the interval 
[fcdec,ni, fcdec,in]- In the vicinity of fcdec.ni, ions stay tightly cou¬ 
pled with neutrals, and their motions are overwhelmed by the 
acoustic perturbations in neutrals. That leads to inefficient 
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(a) WNM 


(b) CNM 



kVA 


^ n 

(c) DC 



U 


(d) SC 


Fig. 1.— vs. kVAll'ni of Alfveti modes in various environments. Solid lines are numerical damping rates of neutral-ion collisional damping. Dashed 

lines are numerical damping rates including both neutral-ion collisional and viscous damping effects. Open circles are analytical damping rates corresponding 
to neutral-ion collisional damping. Triangles show the same as open circles, but for single-fluid approach. The filled circles in (d) show the analytical damping 
rate considering both damping effects. Dash-dotted lines are the cascading rates of Alfven modes. They intersect with the damping rates at fcdam- The analytical 
wavenumbers /cdec.ni. ^dec,in> ^dam> and are denoted by vertical dashed lines. 



Fig. 2.— and \(j^\/uni vs. kVAl'^ni of Alfven modes in MC. 

The same symbols as in Fig. [T]are used, except that the dashed line and open 
squares are the numerical and analytical 


collisional friction and substantially reduced damping rate of 
the ion slow modes. The damping rate is determined by the 
difference between Vnij^ and |a;/| of neutral slow modes (Eq. 
( 15801 ), 

I I Vni Vni . 2 n 2 n 

\uJi\ = — -^ Sin 0 = -^ cos^ 0- (60) 

On the other hand, ions are constrained to magnetic field 
lines. Driven by the magnetic perturbation, the slow modes 
with an effective sound speed 

Cs.eff = \/%Csi = yj2jkBT/mr (61) 

as the phase speed are generated in ions, where the reduced 

mass is = p/rii. It originates from the coupling state 

ions remain with neutrals in this range of wavenumbers. No¬ 
tice that Cs.eff is different from Cs = \/^ksT/rrir of the slow 
modes in strongly coupled regime. The reduced mass in Cg is 
rur = p/n, where n = 2ni + is the total number density. 
It does not distinct the neutral and ionized components since 
neutrals and ions are frozen together and move as one fluid. 
But in the range [A:dec,ni, ^dec.in], ions experience frequent col- 
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(a) WNM (b) CNM 



(c) DC 


Fig. 3.— Same as Fig. [T]but for fast modes. Here we set 9 = 45° as an example. 



kVA 


^ ni 

Fig. 4.— Same as Fig. |2]but for fast modes. 



kVA 


(d) SC 



0 

Fig. 5.— fcdam of fast modes as a function of 6 in WNM. 


lisions with neutrals and hence cannot oscillate freely. Mean¬ 
time, the occasional collisions acting on neutrals are inade¬ 
quate to transfer inertia instantly to neutrals. As a result, we 
see from the expression of Cs,eff that ions alone carry the total 


mass of the medium to move. 

The ion slow modes can only become propagating when 
the real wave frequency = Csefffccos0 becomes larger 
than the damping rate |u;/| (Eq. (l60l l. By applying the cutoff 
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condition |a;i{| = |a;/|, we obtain the cutoff scale 

_ VniCOsO 
^c.tl — 


2c, 


■s,eff 


(62) 


as the lower wavenumber boundary of the ion slow modes 
within [fcdec.ni, fcdec.in]- The ion slow modes emerge at fcc.ti, but 
the phase speed Cs,eff can only be reached above the wavenum¬ 
ber 

Vni cose 

^dec,tl — 7 (63) 

Cs,ef¥ 

satisfying |a;_R| = Vni cos^ 9. It actually signifies ions no 
longer follow the acoustic motions in neutrals, but start to 
develop their own wave modes propagating along magnetic 
field. 

The appearance of the propagating ion slow modes induces 
a new component of the damping rate, associated with the 
wave motions. Since the ion slow modes are driven by mag¬ 
netic perturbation and the propagation is guided by magnetic 
field, no transverse compression can be produced. Therefore 
different from the case for slow modes in strongly coupled 
regime, there is only damping to the parallel propagation. 
Correspondingly, the new component of damping rate is 


Cs,efffe^ COS^ 9 


‘Zi'yi 




(64) 


We now get the complete expressions of the wave frequencies 
of the ion slow modes within [fcdec.ni, fcdec.in] 


= Cs.eff^ cos 6», 

2 n , Cs,etf^^ C0S2 9 


(65a) 

(65b) 


The relative importance of the two terms in Eq. (I65bl i changes 
with k. By equaling the two terms of the damping rate, we 
obtain the transition scale 


h — 

'T'tran — 


Cs.eff 


( 66 ) 


When k < fctran, the first term dominates the damping rate. 
It is independent of uju and corresponds to the decreasing 
\uji\ with k shown in Fig. |6] and Fig. [T] under the condi¬ 
tion of scale-dependent anisotropy. While on wavenumbers 
beyond fctran, the second terms becomes dominant, which is 
proportional to uj\ and increases towards larger k. Notice at 
fctran, |w_r| = Rni COS 9 Can be fulfil led, which coincides with 
the Iw/jI value at fcdec.ni (Eq. (I57al i) and is equivalent to the 
condition Cs,efffc = Rni- We see fctran not only represents the 
transition in damping rate, but also corresponds to the critical 
wavelength within which the effective sound crossing time 
is equal to the neutral-ion collisional time. That is, below 
the scale 1 /fctran, the disturbance associated with the ion slow 
modes cannot be effectively transmitted to neutrals through 
neutral-ion collisions. Neutrals further drift apart from ions, 
and ions are less effectively coupled to neutrals. The differen¬ 
tial motions between ions and neutrals become more signifi¬ 
cant, which consequently causes stronger damping to the ion 
slow modes. 

The propagating component jw/j | starts to decay when 
reaches the value of i/ni, corresponding to the scale 


fcdec,t2 — 


Cs eff cos 9 


(67) 


Neutrals become essentially unaffected by the ion slow 
modes. Subsequently and |a;/| are in equality and the 
ion slow modes are cutoff at the scale 


fcc,t2 — 




Cs.eff cos 9 ’ 


( 68 ) 


which is also the upper wavenumber boundary of the ion slow 

modes within [fcdec.ni, fcdec.in]-^ _ 

The analytical \ujj\ (Eq. (I65bl) ) is shown by open circles 
within [fcdec.ni, fcdec.in] for all the environments in Fig.|6]and|2l 
And [w^given by Eq. (I65al) is also displayed by open squares 
in Fig. |7] The scal es fce ti an d k^^a are indicated by vertical 
dashed lines in Fig. |6(a) |6(d)| The three short vertical dashed 
lines within [fcc.ti, fcc.t2[ in Fig. [T] indicate fcdec.ti, fctran, fcdec.t2 
from large to small scales respectively. Their expressions us¬ 
ing scale-dependent anisotropy are given in AppendixICl 
Table |4] is a summary diagram of the critical scales and 
wave frequencies of ion slow modes, as well as those of the 
slow modes in both strongly and weakly coupled regimes. All 
the expressions apply in low-/3 environment. The scale ex¬ 
pressions present here contain the propagation direction angle 
9. If we disregard the variance of cos 9 with scales, we find 
they are symmetrically linked by 


fctran — 


fcdec.tl 

COS 9 


2k 


'C,tl 


Cnfcc+ 


COS 9 2v^ cos 9 


(69a) 


, . fcc.t 2 Cos 6 » 2^i cos 9k^ 

= fcdcc,t2 COS 9 = , (69b) 


where the last term in Eq. ( I69al ) is obtained by assuming 
sin^ 9^1 and Cs ~ Cgi- If we only focus on the relation 
among the critical scales smaller than k^^^, we find 


2fcdec,t2 = fcc.t2 = -^fcc = -|^fcdec,i 
S71 S71 


(70) 


The same relation also holds for fcdec.ni, k'^, k~, and fcdec.in of 
Alfven and fast modes. That is (see expressions in Table [T] 
and Eq. Q and ( fTbl l) 


2fcdec.ni = kt = ^^fcc = ^T^fcdec.i 


(71) 


when r,,, 1. And again we do not consider the change of 
cos 9 with fc for Alfven waves here. In fact, the wave behavior 
of the ion slow modes below the scale fc,7an fully resembles 
Alfven waves over the whole range of scales we present, with 
VAi replaced by Csi- Notice similar to Cs.eff = y/^Csi, there is 
Ea = V^^VA^. 

Above analysis suggests that with fctran acting as a dividing 
line, the wave spectrum can be divided into two zones, i.e. 
fc < fctran und fc > fctran- The two zones have similarities in 
the sense that the scales fc+, fcc.ti, fcdec.ti in one zone play sim¬ 
ilar roles as fcc.t 2 , k~, fcdec.in in the other. But there also exist 
differences, resulting from the varying coupling state between 
neutrals and ions with scales. As already discussed, one of the 
main differences is for the slow modes in the strongly cou¬ 
pled regime, there is no damping for parallel propagation (Eq. 
(I13bl i'). but for the ion slow modes within [fc^an, fcc,t 2 ], damp¬ 
ing only appears for purely parallel propagation (Eq. 

In addition, we observe |a;/| within [fcc.ti, fctrans] (Eq. (I65bl) ) 
can be approximately written as 


|w/| 


cos^ 9 Vi, 


(72) 
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The factor cos^ reflects much weaker frictional damp¬ 
ing of the ion slow waves compared with jw/l (Eq. dflbl) ) of 
the slow modes in weakly coupled regime. 

To reinforce our understanding on the wave behavior, we 
also provide the force analysis for both fast and slow waves in 
AppendixlP] 

The neutral slow modes originate from the compression 
produced by magnetoacoustic waves. It can only develop af¬ 
ter neutrals decouple from ions and propagate with a sound 
speed in neutral^ In terms of turbulence properties, be¬ 
low the scale l/fcdec,ni, for Alfven and fast modes in low-/3 
medium, hydrodynamic turbulence starts to evolve in neutrals 
where no sound waves arise. But in the case of slow modes, 
neutrals begin to carry acoustic turbulence caused by interact¬ 
ing sound waves (IZakharovll967HZakharov & Sagdeevll970t 
iL’vov et n]|2000l) at scales smaller than l/fcdec.ni- The wave 

motions of neutral fluid also lead to another distinctive feature 
of slow modes from other modes. Instead of having a purely 
nonpropagating interval, propagating ion slow waves appear 
within [/c+j k~]. It suggests that in comparison with hydrody¬ 
namic motions, the wave motions of neutral slow modes cause 
weaker friction against magnetic pressure force on ions, and 
thus the ion slow modes can be driven by the persistent per¬ 
turbations of magnetic held and sustained in ions even at the 
’’cutoff’ wavenumbers. 

Note that in Table|4]and the force analysis of slow waves in 
AppendixiDl we adopt a fixed wave propagation direction and 
disregard its dependence on length scales for simplicity. It is 
necessary to point out, this treatment using fixed propagation 
angle can indeed provide simpler and more intuitively obvi¬ 
ous description of wave behavior than that considering scale- 
dependent anisotropy, but it is only applicable in analysis of 
linear MHD waves. 


The above comparisons with numerical results show that 
the single-fluid approach is able to depict the damping behav¬ 
ior of MHD modes correctly at large scales until reaching the 
cutoff boundary We see for Alfven and slow modes, the 
cutoff boundary /c+ of linear waves is also the fcdam of turbu¬ 
lence in CNM, MC, DC and SC. However, WNM does not 
exhibit a cutoff region for all the three wave modes. The pa¬ 
rameter space required for the existence of wave cutoffs can 
be confined by equaling and k~ . Due to the higher ioniza¬ 
tion degree than other ISM phases, MHD waves in WNM are 
less affected by collisions with neutrals and can avoid being 
cutoff. 

Another remark needs to be made is about the relation be¬ 
tween wave cutoff and fluid decoupling. We found fcdec.ni 
(fcdec.in) and (k~) are of the same order of magnitude, and 
the interval of [fcdec.ni, fcdec.in] is relatively larger than [fc+, k~]. 
Their physical connection is obvious. Only after neutrals sep¬ 
arate themselves from the MHD wave motions of ions, can 
they develop their own motions and exert significant influence 
on ions, i.e. collisional friction strong enough for cutting off 
the MHD waves. On the other hand, although propagating 
waves reemerge at k~, they can only be fully resumed after 
ions get free from coupling with neutrals at fcdec.in- 

6. APPLICATIONS ON CR PROPAGATION IN 
PARTIALLY IONIZED MEDIUM 


^ In high-/? plasma, since the fast modes will behave like sound waves, the 
sound waves in neutrals will become the neutral ’’fast” modes. 


It is necessary to introduce turbulence damping for achiev¬ 
ing a correct unders tanding on the in t eraction of p a rticles and 
MHD turbulence. lYan & Lazariaiil dlOOl 12004 l2008h de¬ 
scribed CR transport in fully ionized medium and clarified 
fast modes are the most effective scatterer of CRs despite their 
damping. It is instructive to extend their study to cover par¬ 
tially ionized medium and build up a complete picture of CR 
propagation in the ISM. 


6.1. Pitch-angle scattering of CRs 

We investigate pitch-angle scattering of CRs based on the 
formalism of di ffusion coeffi c ients obtained from Lokker- 
Planck theory ( lokinii Il966t ISchlickeiser & Milleii 119981: 
lYan & LazarianI 20021 2004) . The energy spectra of dif¬ 
ferent modes we a dopt are obtained from three-dimensiona l 
MHD simulations (ICho et al.l 120021 lYan & LazarianI ^003h . 
Besides, we follow the nonlinear theory (NLT) employed in 
lYan & LazarianI (l2008l) for our calculation on gyroresonance. 
The fluctuations of B in MHD turbulence result in variations 
of particle velocities in both perpendicular and parallel direc¬ 
tions with regards to magnetic held lines. That substantially 
broadens the J-function resonance predicted in quasi-linear 
scattering theory, especially for large-amplitude MHD waves. 
The full set of expressions for Lokker-Planck diffusion coeffi¬ 
cient in low-/3 medium and notations used are presented 
in Appendix lEl (see IVan & Lazarianll2()()4120081 for more de¬ 
tails on the derivation). 

We choose the phases WNM and MC as the representative 
examples of sub- and super-Alfvenic turbulence. Lig. |8]and|9] 
show the calculated 79^^ (normalized by H) as a function of 
fi for 100 PeV CRs in WNM and 30 TeV CRs in MC, where 
/r is cosine of particle pitch angle. CRs’ energy is chosen 
with their Larmor radius larger than the maximum damp¬ 
ing scale among the three turbulence modes. Here 


1 ^. !i = —. 

12 TmpC 


(73) 


We separately calculate the transit-time damping (TTD, 
solid line) for fast and slow modes, and gyroresonance in¬ 
teractions (dashed line) for all turbulent modes. The gyrores- 
onance result calculated using q uasi-linear theory (QLT; see 
iJokipiil [19661 : ISchlickeise^ l200a is also displayed for com¬ 
parison (open circles). Crosses show the approximate result 
using QLT for gyroresonance with fast modes, given by 


77®^ (QLT) 


V'!Ty/JI{l - fi'^) /2 _ 2v^l - 

alVr V7 2V y 


(74) 


This simplified formula is obtained under the condition 
kdam,±TL < 1- We See good consistency between it and the 
numerical integral at small pitch angles. 

The diffusion coefficients in WNM and MC exhibit simi¬ 
lar behavior. Among the three modes, fast turbulence has the 
largest damping scale. Despite severe damping, TTD from 
fast modes dominates CR scattering at large pitch angles. 
On the contrary, gyroresonance is more efficient in scattering 
at small pitch angles. The overall contribution from Alfven 
modes is smaller compared with that from fast modes. Espe¬ 
cially in WNM, due to the more prominent anisotropy over 
all scales in sub-Alfven turbulence, Alfven mo des a re largely 
suppressed in CR scattering. Comparing Lig. |8(b)| and |8(d)[ 
we And fast modes contribute alone in the total diffusion coef- 
fi cient. These results are i n accordanc e with the early findings 
in lYan & LazarianI (l2002Ll^04l2008h . The scattering by slow 
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TABLE 4 

Critical scales and wave frequencies of the slow modes 


Regimes 

Strongly coupled 

ions coupled with neutrals 

weakly coupled 

Scales 

k+ 

^c,tl 

^dec.tl 

^tran 

^dec,t2 

^c,t2 

kc 

^dec,in 

7 COS 9 

Z/t, COS 9 

i/„ i cos 9 

I/tj i 

Uri.i 



Uin 

Cs sin^ 9 

2cs^eff 

Cs,cff 

Cs.eff 

Cg eff COS 9 

Cg eff cos 9 

2csi cos 9 

Csi cos 9 

I'^hI 

< Csk COS 6 

< Cs,efffc COS 9 

Cs,eff^ COS 9 

< Cs,eff^ COS 9 

< Cgik COS 9 

Csik COS 9 


sin^ e 

j cos^ 9 1 *^s,eff^ ^ 

2 * 

Vin 



(a) WNM 



(c) DC 



(b) CNM 



(d) SC 


Fig. 6. — Same as Fig. [T]but for slow modes. Only one branch of slow modes is present in strongly coupled regime {k < fcdec,ni)- On wavenumbers beyond 
^dec,ni 5 there are both neutral and ion slow modes. Their analytical damping rates are represented by filled and open circles respectively. Four cutoff boundaries 
> ^c,£i 5 ^c ,£2 indicated in the figures. 


modes depends on /3 of the medium (see Eq. (IE16l) - (IE181 l). 
Low /3 value in ISM conditions (see Tableful makes the effect 
of slow modes negligible. In addition, by comparing with the 
results of gyroresonance from QLT, we find nonlinear effect 
results in a much wider range of pitch angles, while QLT gives 
D/Mfi = 0 at large pitch angles due to its discrete resonances, 
fc||,res = This resonance condition cannot be satisfied 

at small p, as the turbulence is damped before cascading down 
to small scales. Only at the large p end, gyroresonance from 
QLT comes into coincidence with that from NET. 


To examine the damping effect on CR scattering, we again 
calculate the diffusion coefficients but for CRs with relatively 
low energies. CRs’ energy is chosen with their Larmor radius 
rL smaller than the minimum damping scale among the three 
turbulence modes. But since in WNM, slo w m odes survive 
neutral-ion collisional damping (see Section ISJT l and can ex¬ 
tend to the gyroscale of ions, we choose 10 TeV CRs with 
smaller than the damping scale of Alfven modes, which is 
smaller than that of fast modes. Eor MC, Alfven modes have 
the smallest damping scale. Accordingly, we adopt 10 GeV 
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V 


Fig. 7.— Same as Fig. |2]but for slow modes. The filled symbols con'espond 
to the neutral slow modes. The shoif vertical dashed lines within [fc^.n, ^ 0 , 12 ] 
represent fcdec.tl. ^tran, fcdec,t2 from left to right respectively. 

as the CR energy. 

Fig. [To] and Fig. [TT] show the results in WNM and MC 
respectively. The asterisks represent the simplified expression 
for TTD with fast modes, 

^ exp 

(75) 

which is obtained under the condition fcdam,_LTL < 1 and 
agrees well with integral result for low-rigidity CRs over all 
pitch angles. Obviously, TTD with fast modes becomes the 
only significant scattering effect. The remaining gyroreso- 
nance of slow modes in WNM has negligible effect. Com¬ 
pared with the results for high-energy CRs with larger than 
the damping scales of turbulence modes in Fig. |8|and|9| here 
we can clearly see the resonance gap at small pitch angles due 
to the absence of gyroresonance. TTD alone cannot fulfill the 
scattering of CRs with small pitch angles. 


6.2. Scattering frequency of TTD and gyroresonance 

To gain a deeper understanding on the dependence of TTD 
and gyroresonance on CR energy and the influence of turbu¬ 
lence damping, we illustrate the scattering frequency v = 
— fi^) of TTD and gyroresonance separately in 
WNM (Fig.fTZt and MC (Fig. [T^ . Since the major scattering 
agent differs at small and large p itch angles , we present results 
of TTD at /r = 0.1 in Fig . |12(a)|and|13(a)| and gyroresonance 
at p. = 0.8 in Fig. 1 12(b) and 1 13(b)] The vertical dashed lines 
indicate the CR energies with equal to the damping scales 
of Alfven (i?dam, a) and fast (i?dam, f) modes. Slow modes are 
not considered due to their n egligib le effect. 

In WNM, as shown in Fig. |12(a) when damping is absent, 
1 / (dashed line) through TTD with fast modes decreases with 
CR energy. Otherwise 1 / (solid line) keeps steady until reach¬ 
ing i?dam, f (vertical dashed line), and then coincides with the 
result without damping. The analytical scattering frequency 
(asterisks) is derived by using Eq. ( ITST l, and can apply to 
the energy range below i?dam, f- As illustrated by Eq. (fTSl l, 
for CRs with is determined by the damping 

scale. Since CR velocity is approximately equal to the light 
speed and its change with energy is marginal, Eq. ( l75l l yields 
a constant and i/ at a fixed p. It implies only the turbu¬ 


lence on scales larger than contributes to TTD interaction. 
As a result, in the case of no turbulence damping, 1 / of lower- 
energy CRs has a higher value because turbulence on a larger 
range of scales is involved in TTD scattering. While in the 
case with damping, ly of TTD is independent of CR energy 
for CRs with < l/fcdam, and decreases with energy when 

tl > 1/fcd am- 

Eig. |12(b)| shows the total 1 / of gyroresonance by both 
Alfven and fast modes in WNM. The results free of damping 
and in the presence of damping come into coincidence when 
CR energy is over E’dam, f- The respective roles of Alfven and 
fast modes are explicitly displayed. When there is no damp¬ 
ing, Alfven modes have negligible effect over the whole en¬ 
ergy range due to strong turbulence anisotropy in WNM. So 
o for gyroresonance with fast modes (open triangles) coin¬ 
cides with the total ly. The crosses are analytical ly by em¬ 
ploying Eq. (TTTI i. Although Eq. (l74l i is the simplified 
using QLT for fcdam,_LfL < 1, the gyroresonance condition 

k\\ = n/(up) = ^/T — p^/(rip) can still be satisfied at 
large p, so that it can reach a good agreement with the nu¬ 
merical result using NLTQ Similar to TTD, gyroresonance 
with fast modes has decreasing ly with increasing CR ener¬ 
gies, but their physical origins are different. Eor TTD, there 
is not a specific resonant scale. Turbulence over all the scales 
above tl contributes to particle scattering. However, gyrores¬ 
onance can only become effective on specified scales which 
are determined by CR energy. We observe in Eq. (i74l) that 
oc IjVR, where R is CR rigidity. It demonstrates 
that higher-energy CRs are scattered by larger-size turbulence 
eddies, and consequently have larger mean free paths and 
smaller and ly. 

When damping is taken into account, since gyroresonance 
has a preferential resonant scale, 0 corresponding to Alfven 
(filled circles) and fast (filled triangles) modes only become 
efficient after exceeds their damping scales. The analyt¬ 
ical given by Eq. (l74l l can only apply to CRs with en¬ 
ergies larger than i?dam,f in this situation. Therefore, TTD is 
the only effective interaction and governs particle scattering at 
lower energies. Since Alfven modes have a smaller damping 
scale than fast modes, their contribution can be seen within 
the energy range [i?dam,A) Fidam.f]- But fast modes dominate 
the scattering of CR with energies higher than Eidam.f- 

In the case of super-Alfvenic MC, we see much higher v 
values for both TTD and gyroresonance. The ly of T TD in Eig. 
|13(a)| shows a similar trend as tha t in Ei g. |12(a) for WNM. 
But for V of gyroresonance (Eig. |I3(b)| l, the Alfven modes 
become more important in MC compared w ith WN M due to 
less prominent turbulence anisotropy. Eig. |13(c)| and |13(d)| 
display the components of v from Alfven and fast contribu¬ 
tions for cases without and with damping. When damping is 
absent, ly of gyroresonance with Alfven modes (open circles) 
increases towards higher energies. The opposite trend com¬ 
pared with gyroresonance with fast modes originates from the 
scale-dependent turbulence anisotropy. Lower-energy CRs in¬ 
teract with many elongated eddies with k± ^ A:|| ^ l/r^ 
within one gyroperiod and this random walk causes inefficient 
scattering. While higher-energy CRs are scattered by larger- 
size eddies which are less anisotropic and he nce the scattering 
efficie ncy increases (lYan & Lazarianll2004ft . We see in Eig. 
|13(c)[ approaching Za, weak turbulence anisotropy makes 

^ In fact, Eq. llU is only valid for > 0.5314 to ensure a positive value 
of Dan. 
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Fig. 8 .— Pitch-angle diffusion coefficients (nonnalized by O) of 100 PeV CRs in (a) Alfven, (b) fast, and (c) slow modes in WNM. Solid and dashed lines refer 
to TTD and gyroresonance hy applying NLT. Open circles are QLT results for gyroresonance. Crosses in (b) show the approximate QLT result for gyroresonance 
given by Eq. 1741 . (d) Total diffusion coefficients from three modes. The dash-dotted line is the sum of the solid and dashed lines. 


Alfven modes have comparable efficiency as fast modes in 
gyroresonance scattering. Their joint effect flattens v at high 
energies. 

In brief, we And both TTD and gyroresonance with fast 
modes decrease with CR energy, while gyroresonance with 
Alfven modes increases with energy. The damping effect can 
only make a difference in CR scattering when < l/Aidam- 
Then TTD becomes the only scattering agent and has a scat¬ 
tering efficiency independent of particle energy. At higher 
energies with > l/Aidam, the importance of gyroresonance 
with Alfven modes depends on the anisotropy degree of the 
turbulence. In super-Alfvenic turbulence, Alfven modes can 
become more efficient than fast modes in gyroresonance scat¬ 
tering at large scales comparable to Ia- 


6.3. Results on parallel mean free path of CRs 

Given the results of the parallel mean free path Ay of 
CRs can be evaluated by integrating over p,. 


3v 

T 



(I-mT 


(76) 


where is the total contributions from both TTD and 
gyr oresonance w ith the three modes calculated using NLT. 
Fig. 1 14(a^l4(d)| present Ay as a function of CR energy Ek- 


The upper limit of the energy range is restrained by < L 
in WNM and < lAtn other environments. 

We can clearly see the imp act of turbulence damping on 
Ay. In the case of WNM (Fig. |14(a)| l, even in the undamped 
range of turbulence, neither Alfven nor slow modes can effec¬ 
tively scatter CRs. Only when CR energy approaches i?dam, f, 
can CRs be confined by gyroresonance with fast modes. Nev¬ 
ertheless, the resulting Ay is larger than Itr over the whole 
energy range. It shows in WNM, direct interactions with tur¬ 
bulence modes are invalid in scattering CRs. Instead, other 
mechanisms, e.g. held line wandering, can play a more im¬ 
portant role in determining Ay. _ 

In the other ISM phases (Fig. |14(b)fT4(d)l l, we observe 
a nonmonotonic U-shaped dependence of Ay on CR energy. 
The dramatical decrease of Ay at ii^dam, a and the further de¬ 
crease of Ay at i?dam, f come from the contribution of gyroreso¬ 
nance with Alfven and fast modes respectively, while the rise 
of Ay at Ek > i?dam, f is in accordance with the decreasing 
1 / of TTD with energy. The marginal change of Ay at E'dam, s 
reflects the insignificant role of slow modes in CR scattering. 

On the other hand, although TTD contributes over all en¬ 
ergy range and is more efficient in CR scattering, it is inca¬ 
pable to scatter CRs with small pitch angles. Thus TTD alone 
is not enough to confine CRs and leads to infinite Ay in the 
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Fig. 9.— Same as Fig.[^but for 30 TeV CRs in MC. 




Fig. 10.— Same as Fig.[8]but for 10 TeV CRs in (a) fast and (b) slow modes in WNM. The asterisks in (a) are analytical result given by Eq. (75). 
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Fig. 12.— Scattering frequency vs energy of CRs in WNM from (a) TTD at /i = 0.1 and (b) gyroresonance at /j. = 0.8. Dashed and solid lines correspond 
to cases in the absence of damping and in the presence of damping. Circles and triangles in (b) show the separate contributions from Alfven an d fa st mo des. 
Filled and open symbols correspond to situations with and without damping. Asterisks and crosses represent the analytical results by using Eq. and o 
respectively. 
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(c) Gyroresonance, without damping 



Fig. 13.— Same as Fig. ll2l but for MC. (b) is replotted in (c) and (d), showing the separate contributions from Alfven (circles) and fast (triangles) modes. 


energy range where gyroresonance is absent. Therefore, the 
lower energy limit of effective CR scattering and finite Ay is 
determined by the damping scale of Alfven mode. 

Table |5] summarizes the lowest energies of CRs with their 
scattering unaffected by turbulence damping in different par¬ 
tially ionized ISM, which are determined by the damping 
scales of fast modes. The scattering and acceleration of 
CRs with lower energies are subject to turbulence damp- 
ing. Notice that diffe rent from the situation discussed here, 
lYan & Lazarianll2004l showed in fully ionized plasma, damp¬ 
ing of fast modes strongly depends on the angle 9 between k 
and B, which results in anisotropic distribution of fast mode 
energy at small scales. In their case, the lowest energy for 
CR scattering unaffected by turbulence damping should be a 
function of 9. 

6.4. Other mechanisms on confining CRs 

Our results show that in partially ionized medium, MHD 
turbulence, especially fast modes, is severely damped by 
neutral-ion collisions. Consequently, only high-energy CRs 
can be efficiently scattered through interactions with turbu¬ 
lence modes. Besides scattering, the diffusion of CRs can 
also arise from the spatial wandering of magnetic field lines 
(lJokir)iilll966l) . In three dimensional turbulence, the wander- 


TABLE 5 

The minimum energy oe CRs unaeeected by 

TURBULENCE DAMPING IN DIEEERENT ISM PHASES. 


ISM phases 

^dam 

^k,min 

Alfven 

fast 

slow 

WNM 

0.003 pc 

4.0 pc 

— 

45.3 PeV 

CNM 

0.005 pc 

0.1 pc 

0.04 pc 

1.2 PeV 

MC 

6.7 AU 

0.002 pc 

98.2 AU 

18.9 TeV 

DC 

35.0 AU 

0.009 pc 

261.7 AU 

0.99 PeV 


ing of field lines is induced by t urbulent motions and charac- 
terized by turbulence properties. iLazarian & Vishniak (119991) 
quantitively described the field line wandering according to 
the scaling laws of GS95-type MHD turbulence, and demon¬ 
strated its key role in determining th e rate that magnetic 
reconn ection proceeds. Based on the ILazarian & VishniacI 
(Il999h prescription for magnetic field wandering and the re¬ 
lated field line diffusion, signi ficant theoretical reformula- 
tions on e.g. thermal production (iNaravan & Medvedevl200i] : 
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Fig. 14.— Parallel mean free path of CRs as a function of their energies in (a) WNM, (b) CNM, (c) MC and (d) DC. Vertical dashed lines represent the CR 
energies with their equal to the damping scales of Alfven (F^dam, a )5 slow (F^dam, 5)5 fast (F^dam, f) modes. Notice that F^dam, A is below the energy range 
shown in (a). Etr is the energy con'esponding to Itr- The horizontal dashed lines in (b), (c) and (d) refer to the length scales of 


Lazarianll2006l) and CR pro pagation (lYan & Lazarianll2008L 
2Q12HLazarian & Yanll2014h . have been achieved. 

Other mechanisms which can enhance sca t tering of CRs in- 
clude streaming instability (IWentzellll974t ICesarskv 198C) 
and gyroresonance in stability (iLazarian & Beresnvak 2006 : 
lYan & Lazarianll201 Ih . The additional MHD waves excited 
by CRs can in turn efficien tly scatter CRs w ithin a range of 
energies (e.g., < 100 GeVI. lYan et al.l (1201 2h investigated the 
shock acceleration of CRs in the presence of streaming in¬ 
stability. They found the CR flux near supernova remnants 
is strongly enhanced compared with typical Galactic values, 
so the growth rates of streaming instability can overcome the 
background turbulence damping and boost scattering of CRs. 
This finding enables them to reproduce the observed gamma- 
ray emission from the supernova remnant W28. 

The growth rate of the resonant waves excited by streaming 
CRs is 


r&cR 


ncR(> r) 

iZo- 

m 


^Stream 

Va 



(77) 


where Hq = eB/mpC is the non-relativistic CR gyrofre- 
quency. ncR(> T) is the number density of CRs with their 
energies larger than T GeV, and the value correspond¬ 
ing to r gives the wavevector of the generated waves, i.e. 
^ l/fc|| 1/fc, which propagate closely parallel to mag¬ 


netic field. And ^stream is the streaming velocity of CRs, which 
should be larger than Va so as to amplify waves. 

The instability can only set in when Tg.cR exceeds the to¬ 
tal damping rate, which includes the neutral-ion collisional 
and neutral viscous damping discussed in this paper, and 


also the nonlinear dami 

ping by background MHD turbulence 

(lYan & LazarianI 12002 

12004 

t iLarmer & Goldreichl 12004 

IBeresnvak & LazarianI 

200i" 

Yan & LazarianI 2011). The 


non-l inear damping rate takes the form as (lYan & LazarianI 

EMI 


T,urb=rL^lA^VA (78) 

for tl < Za in super-Alfvenic turbulence, and 

Tturb =rlh-iVL (79) 

for < hr in sub-Alfvenic turbulence. 

As an example. Fig. [15] compares the growth and damp¬ 
ing rates at different CR energies in the WNM environment. 
The neutral-ion collisional damping rate for parallel Alfven 
waves (Eq. ( l4b] i with 0 = 0) applies, which has the value 
of jw/l = r'm/2 at large wavenumbers corresponding to low 
CR energies. Notice neutral viscous damping is unimportant 
for the generated waves in this particular case. The nonlinear 
damping only exists for CRs with their larger than the par- 
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allel damping scale of the background Alfvenic turbulence. 
For the growth rate represented by the solid line, we assume 
t^stream has the Same order of magnitude as Va and adopt the 
CR density ne^ the sun, ncR(> F) = 2 x 10“^°r“^ ®cm“^ 
(lWentzellll974ft . It shows the streaming instability can only 
overcome the collisional damping and effectively contribute 
to particle scattering for CRs with energies lower than ^10 
GeV. Crosses show a higher growth rate with the CR density 
increased by three orders of magnitude. As a result, CRs with 
up to ^ 1 TeV can be significantly scattered and conhned. 



Fig. 15. — Growth and damping rates (normalized by i/„i) of streaming in¬ 
stability as a function of CR energy in WNM. The dashed line is the neutral- 
ion collisional damping rate. The dash-dotted line is the nonlinear damping 
rate. The growth rates with different CR number densities, ncR(> F) = 
2 X 10“^''r“^'®cm~® (solid line), and ncR(> F) = 2 X 10“^F“^'®cm“® 
(crosses) are shown. The vertical dashed line indicates the energy come- 
sponding to the parallel damping scale of the background Alfven modes. 

From above example we see to establish a comprehen¬ 
sive picture of CR propagation, other scattering mechanisms 
should also be incorporated in addition to direct interac¬ 
tions with turbulence modes. For relatively low-energy CRs, 
streaming instability can be a promising mechanism of CR 
confinement by scattering. For high-energy CRs with their 
exceeding the damping scale of MHD turbulence, TTD and 
gyroresonance scattering operate. Regarding the CRs with in¬ 
termediate energies, their confinement is mainly attributed to 
field line wandering as mentioned earlier. 

We again stress that in analyzing these scattering processes 
in a partially ionized medium, neutral-ion collisional damping 
is an essential physical ingredient, which determines the dissi¬ 
pation scale of turbulence cascade, and is the major constraint 
for the growth of instabilities. Taking the damping effect into 
account is necessary to attain realistic calculation of diffusion 
coefficients and proper understanding of CR propagation in 
partially ionized medium. 

7. DISCUSSION 

In partially ionized medium, the wave behavior strongly 
depends on the coupling state between neutrals and ions. 
We consider low-/3 medium and specifically distinguish the 
neutral-ion and ion-neutral decoupling scales of three modes 
and the regimes with different coupling degrees separated by 
them. The cutoff boundaries of waves are closely correlated 


with the decoupling scales. In fact, the cutoff interval approxi¬ 
mately shares the same domain as that between the neutral-ion 
and ion-neutral decoupling scales. This coincidence indicates 
fluid decoupling can be the physical origin of the wave cutoff. 
Unlike Alfven and fast waves, propagating waves can arise 
within the original cutoff region of slow waves, which results 
in two cutoff intervals sitting among three propagating wave 
branches. We for the first time provide full analytical expres¬ 
sions of all the cutoff boundaries and the wave frequencies 
of both ’’neutral” and ”ion” slow waves in different coupling 
regimes. Slow waves on scales within [fcdec.ni) fcdec.in] exhibit 
more complicated properties in both propagating and damp¬ 
ing components of the wave frequencies. Accordingly, the 
wave spectrum on intermediate scales is also divided into dif¬ 
ferent coupling regimes, conhned by multiple critical scales 
whi ch are linked together in a symmetric patten (see Section 

O. 

In order to test the validity of the single-Huid approach, we 
adopted Newtons iteration method to solve the single-huid 
dispersion relation. By comparing the solutions with those 
to the two-huid dispersion relation under weak damping as¬ 
sumption, we see consistent results at large scales down to 
the lower cutoff boundary We further numerically con- 
hrmed the validity of single-huid approximation in describing 
MHD waves in partially ionized ISM and SC environments. It 
means for practical purposes, single-huid treatment can also 
be used to determine the damping scales of both incompress¬ 
ible and compressible turbulence modes in strongly coupled 
regime in similar conditions. 

To obtain the damping scales of MHD turbulence, we fol¬ 
low the present-day understanding of turbulence and treat the 
turbulence cascades of Alfven, fast and slow modes sepa¬ 
rately. It turns out the scale-dependent anisotropy of turbu¬ 
lence plays a critical role in regulating wave behavior and de¬ 
riving turbulence damping. We see a close relation between 
the cutoffs appearing in linear MHD waves and the damping 
of MHD turbulence. Particularly, because of the intrinsic crit¬ 
ical balance of the turbulent motions (GS95), the cascade of 
Alfven modes is truncated at the lower cutoff boundary where 
the nonpropagating waves arise. As a result, the wave cutoff 
should also be taken into account when deriving turbulence 
damping scales. 

In addition, we apply the analytical results on turbulence 
damping to a variety of partially ionized ISM phases and solar 
chromosphere. 

- We hnd neutral viscosity plays a significant role in 
damping Alfven modes in WNM and SC, while neutral- 
ion collisions act as the dominant damping effect in 
other environments. 

- Fast modes in all conditions are damped out in strongly 
coupled regime. And the damping is especially severe 
in WNM and SC because the cascading rates in both 
environments are substantially low and in addition, the 
damping rate in WNM is relatively high due to its low 
density. Besides, different from the case in fully ionized 
medium, the damping of fast modes only has a weak 
dependence on wave propagation direction. 

- In the case of slow modes, since cutoff appears earlier 
than turbulence damping, the damping scales of slow 
modes are given by the largest cutoff scales in all condi¬ 
tions except for WNM. Due to the absence of cutoff in¬ 
tervals, slow modes in WNM survive neutral-ion colli- 
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sional damping. In fact, all the three modes in WNM do 
not have cutoffs. The physical explanation is due to the 
relatively high ionization fraction in WNM, coupling 
of neutrals with ions on larger scales is easier since ev¬ 
ery neutral has more chances to collide with ions, and 
decoupling of ions from neutrals on smaller scales is 
also easier since every ion has fewer chances to collide 
with neutrals. Consequently, the separation between the 
neutral-ion and ion-neutral decoupling scales is much 
shortened and cutoff can be avoided. 

Regarding the turbulence cascade of slow modes below the 
neutral-ion decoupling scale in WNM, although the MHD 
turbulence in ions survives neutral-ion collisional damping, 
since slow modes are slaved to Alfven, it may be quenched at 
a smaller scale with the damping of Alfven modes. Provided 
the cascade of ”ion” slow modes can extend to the regime 
w here collisionless damp i ng is d ominant, we refer the reader 
to iLithwick & Goldreichl (1200 Ih on the slow mode damping 
below the proton mean free path. For the acoustic turbu¬ 
lence of’’neutral” slow modes, it is subject to both neutral-ion 
collisional and neutral viscous damping. The damping scale 
can be determined by comparing the cascading rate and both 
damping rates, which falls beyond the scope of the current 
work and will be discussed in future work. 

As one important application of the present study, we ex¬ 
plored the damping effect on CR propagation in partially ion¬ 
ized medium. The local reference system where GS95-type 
scalings stand is also very important for studying CR propa¬ 
gation, since the scattering of a CR particle is determined by 
its interaction with the local magnetic field perturbation in¬ 
stead of background held. That ensures realistic statistics of 
magnetic huctuations which we use for calculating CR scat¬ 
tering. 

Al fven modes are f ound to be an ineffec t ive sc atterer of 
CRs (IChandranll2000t lYan & Lazarianll2002L l2004t) . During 
the energy cascade of Alfven modes, turbulent eddies become 
more and more elongated along magnetic held lines, with en¬ 
ergy concentrating in the direction perpendicular to the local 
mean magnetic held. Thus a CR particle with compara¬ 
ble to the parallel scale of eddies interacts with many uncor¬ 
related eddies in perpendicular direction within one gyrope- 
riod. The random walk leads to very inefficient gyroresonance 
scattering of CRs by Alfven modes. But in the presence of 
neutral-ion collisional damping, turbulence cascade is heav¬ 
ily damped and truncated at a large scale, where turbulence 
anisotropy is relatively weak. Especially in super-Alfvenic 
turbulence, turbulent eddies are nearly isotropic when ap¬ 
proaches I A- As a result, high-energy CRs can be effectively 
scattered by gyroresonance with Alfven modes, as shown 
above. 


8. SUMMARY 

In this paper we continued the work in Paper I on damp¬ 
ing of MHD turbulence in partially ionized medium. We put 
more emphasis on compressible modes and apply the analyti¬ 
cal results to a variety of low-/3 ISM phases and solar chromo¬ 
sphere. As an important application, we studied the CR prop¬ 
agation in partially ionized medium considering the neutral- 


ion collisional damping, as well as neutral viscous damping, 
of MHD turbulence. Here we summarize the main results. 

1. We present explicit correlation between the cutoff 
boundaries of MHD waves and decoupling scales. Under the 
consideration of scale-dependent anisotropy, they are 


kt 


k 


C 


for Alfven modes, and 


( > ) ^ ^dec.ni: 

2 ^ ^dec.in- 


kt 

kc 


^ ^dec.ni; 
Sn 

i*dec,in. 


for fast modes. Slow modes have the same 

kf, =2 ^ Ajdec.in 


(80a) 

(80b) 

(81a) 

(81b) 

(82) 


as Alfven modes, but the relation between k'^ and fcdec.ni de¬ 
pends on turbulence properties. 

2. Single-fluid approach is capable of correctly describ¬ 
ing wave behavior and damping properties in strong coupling 
regime. 

3. Cutoff of MHD waves should also be taken into account 
when deriving damping scales of turbulence cascade. 

4. We showed the importance of neutral viscosity in damp¬ 
ing Alfven modes in WNM and SC. Especially for the Alfven 
modes in SC, neutral viscosity is the dominant damping effect 
instead of neutral-ion collisions. 

5. We thoroughly investigated the behavior of slow waves 
and provided full expressions of the multiple cutoffs and wave 
frequencies in different coupling regimes. 

6. The scale-dependent anisotropy of GS95-type turbulence 
is important for understanding wave behavior, calculating tur¬ 
bulence damping scales, and studying interactions between 
CRs and magnetic perturbations. 

7. We evaluated the scattering efficiencies of TTD and gy¬ 
roresonance with the three turbulence modes, and found due 
to the severe damping in partially ionized medium, only high- 
energy CRs can be effectively scattered through direct inter¬ 
actions with turbulence modes. Compared with other environ¬ 
ments, CRs in WNM are poorly confined due to its prominent 
turbulence anisotropy, as well as the largest damping scale of 
fast modes. As for CRs with lower energies, other effects such 
as held line wandering and streaming instability can set in and 
contribute in confining CRs. 
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drag coetticient 

Id 

ion density 

pi 

neutral density 

pn 

total density 

p 

proton mass 

rup 

mass of hydrogen atom 

mu 

ion mass 

mi 

neutral mass 

rrin 

reduced mass 

mr 

ion number density 

Ui 

neutral number density 

rin 

total number density 

n 

neutral-ion collisional cross-section 


electronic charge 

e 

neutral-ion collision frequency 

^ni 

ion-neutral collision frequency 

^in 

ion fraction 


neutral fraction 


ratio of neutral and ion densities 

X 

wave frequency 

LO 

real part of wave frequency 


imaginary part of wave frequency 

LJI 

wave number 

k 

x^y, z component of k 

kx ■jky^kz 

k component parallel to local magnetic field 

h 

k component perpendicular to local magnetic field 

fcx 

wave propagation angle with regard of magnetic field 

9 

Alfven speed 

Va 

Alfven speed of ion-electron gas 

VAi 

magnetic field 

B 

sound speed 

Cs 

sound speed of neutral gas 

C-sn 

sound speed of ion-electron gas 

^si 

ion gas pressure 

Pi 

neutral gas pressure 

Pn 

temperature 

T 

Boltzmann constant 

ks 

adiabatic constant 

7 

ratio of gas pressure to magnetic pressure 

P 

injection scale of turbulence 

L 

injection scale of MHD turbulence in super-Alfvenic turbulence 

Ia 

transition scale from weak to strong MHD turbulence in sub-Alfvenic turbulence 

hr 

turbulent velocity at L 

Vl 

turbulent velocity at 1 

VI 

turbulent velocity at hr 

Vtr 

Alfvenic Mach number 

Ma 

damping scale 

Id 

ion velocity 

Vi 

neutral velocity 

Vn 

mean free path for a neutral particle 


cross section for a neutral-neutral collision 

O' nn 

speed of light 

C 

Lorentz factor 

r 

CR particle’s velocity 

V 

CR pai'ticle’s gyrofrequency 

Q 

CR pai'ticle’s non-relativistic gyrofrequency 

f!o 

CR particle’s pitch angle cosine 

P 

pitch angle diffusion coefficient 


scattering frequency 

u 

growth rate of streaming instability 

rg,cR 

streaming speed of CRs 

'^^stream 

number density of CRs 

ncR 

nonlinear damping rate 

r turb 


APPROXIMATE SOLUTIONS TO SINGLE-FLUID DISPERSION RELATION 


Here we list th e ana lytical solutions to the single-fluid dispersion relation at some tractable limits following the method de¬ 
scribed in Section lU^ 

1 . 0 ° 


2.0^ 90° 


Ai ,2 = ±c„,A3,4 = —{-ik± y4- k^). 


Ai ,2 = ±c„. 


(Bl) 

(B2) 


3. /3-^ 0 


Ai, 2 = ±Cn COsO 

A 3,4 = X^{—iksec9 ± \/4 — sec^ 9). 


(B3) 
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4 . /3 —>■ oo 

iVACOs6{ksec^ 9 ± \/k^sec‘^0 — 4) 
Ai, 2 = ±c„, A3,4 =--- 


5. fc —>■ 0 (ideal MHD case) 


where 

6. Small 0 and /3 —0 


where 

7. Small 0 and /3 —>■ 00 


8. 6» < 90 ° 


9. Low /? and 0 —> 0 

10. Low /3 and 0 ^ 90° 

11. High/3 and 0 ^ 0 

12. High /3 and 0 ^ 90° 

13. Small k and /3 —>■ 0 

14. Small k and /3 —>■ 00 


A 1.2 = ± 


y'cl + Fj + A 

V2 


A 3,4 = ± 


\/ Cp- 
V '-'n 


yj-A 


72 


A = ^ (4 + 71)^ - 44yj cos2 0 . 

Ai,2 = 

A, = -lAC,- _ 4C-2fcsecec"+C^ „■ 

2 4(2C2+3fcsec^C-4) 

A • 4k— 2fc sec A 

^ 2 4(2^2 —3A: sec 4) 


C = k — iy 4, — kP-^n = k-\-iy A — kP-. 


Ai.2 = ±C„ 
A3 = -^* 
A 4 = 


4 COS^ ^C + C^ .7 

4(fc sec ^ —C) 

4 cos^ 9k.^k^ 

4{k sec 9 — k) 


— Cn 


Cn sin^ 9Vi 


2(c^+ZCnfc sec OVA — yX) ^ 


A 2 — Cn 


A3 = — 


fc sec^ 0 


2(c^ —fcnfc sec 9Va — VX ^ 


. , . VaC . , Vak. 

Ai,2 = ±C„, A 3 =- —I, A 4 =- 


Ai, 2 = -^(—ifctan0 ± \j A — k^ tan02). 


, _ y4(fc±7P-4), 

^1,2 — ^^72, A3^4 — 2. 


Ai. 2 = ±c„, A3 = —iVAk sec 0. 


7,2 =±^, 

\ _ TT- ikVA sec 9 \ _ T/ i feVAsec^ 

A3 — — 3 ^fcsece+ 2 ’ ^ ^ 3 fcsec 6 »+ 22 ' 


A1j2 — ^Cn-, 

\ _ Va ikVA sec 9 \ _ Va 

'^3 — /?: r^/ •, /s , /7:\ 5 ^4 — /tt 


a/2 \pl{ik sec 0 +a/2) ' 


ikVA sec 9 


vA y2( — f fe sec e+y2) ■ 


15. Small k and 0 —>■ 0 


A 3 


Va + 


tkVl(Vl-cl) 

2kVA{cl-3V^)+2VA{cl 


Ai,2 — 

, A 4 


-Vi) 


—I—^ 

= —Va + 


Vkvi{vi-cl) 

-ikVA{cl-3Vl)+2VA{cl-Vl) ■ 


16. Small k and 0 —90° 


Ai,2 


^ 2(yj + c^)i 

371 + 24 


(B4) 

(B5) 

(B6) 

(B7) 

(B8) 

(B9) 

(BIO) 

(BID 

(B12) 

(B13) 

(B14) 

(B15) 

(B16) 

(B17) 


(B18) 
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CRITICAL SCALES OF ION SLOW MODES CONSIDERING SCALE-DEPENDENT ANISOTROPY 
By substituting the scaling relations described in Eq. (l30l l and (l32l l into Eq. (l6^ . db^ . dCTl i. and dbSl) . we get expressions 


fccti = f , ^ ^h,2 = 


^dec,tl 


V ) 


2z/r7 




V 

for 1/fc < IA in super-Alfvenic turbulence, and 

/ 

/Cftl = I - 


L 1 ^dec,t2 — 


^dec,t 


\2V5’c 

_ f ^ni \ 

V y/^iCsi ) 


L ^Ma, kc,t2 


= (- 


's/^Csi 

^ni 


vVS’c 


L ^ Ma , ^dec,t2 — 


VTic 


3 

2 1 -3 






for 1/k < Itr in sub-Alfvenic turbulence, under the consideration of strong turbulence anisotropy. 


(Cl) 

(C2) 

(C3) 

(C4) 


FORCE ANALYSIS OF MAGNETOACOUSTIC WAVES 

The wave behavior can be better understood with the aid of force analysis. Here we follow the method in ISoler et al.l (12013bl) 
for Alfven waves to perf orm an analysis on the fo rces for magnetoacoustic waves. We start from the momentum equations of 
ions and neutrals (see e.g. iZaaarashvili et al.ll20Hl) . 


-ioJVix = 

^2 iL2 
A^si^x 1 

[Vix +Viz) ■ 

:‘^Vlkl 

'^ix 

UJ 

UJ 


^inip^ix '^nx)-) 



—iojViz = 

^2 iL2 

^ ^si^x j 

UJ 

P^ix '^iz') 

~ ^iniyiz ~ 

'^nz)-) 

iujVfix — 

„2 l2 
^sn^x 

i'^nx ‘^nz 

) “1” ^nii.'^ix 

"^nx 

UJ 

-iujVnz = 

: ^sn^x 

i'^nx '^nz. 

) ^nii^iz 

- Vnz[ 

UJ 


(Dla) 

(Dlb) 

(Die) 

(Did) 


Here Vi and are velocity perturbations in ions and neutrals. We consider the constant magnetic field is along z-direction, and 
waves propagate in x-z plane. For a qualitative study on wave behavior, we fix the wave propagation direction as 45° with regards 
to magnetic field, i.e. kx = kz, to simplify the algebraic procedures. To conduct the force analysis for ion and neutral particles in 
different directions, from above equations we get 


= |P,| + |M|-|F,|, 
fi,z = \PX - \Fz\, 
fn,x = \Pn \ + \Px\, 
fn,z = \Pn\ + \Fz\. 


(D2a) 

(D2b) 

(D2c) 

(D2d) 


Here \Pi\ and |P„| are gas pressure forces acting on ions and neutrals. \M\ is magnetic pressure force and only acts on ions in 
perpendicular direction to magnetic field. \Fx \ and \Fz \ are friction forces between neutrals and ions in x and z directions. The 
expressions of forces are not displayed for simplicity. We next take the parameters of MC phase and use the numerically solved 
oj of fast and slow waves to calculate the corresponding forces. To further simply the problem, we also set rrii = 2m„ to have 
equal sound speeds in neutrals and ions. 

In figure[T6] we display the wave frequencies and different forces |/| normalized by \M\ for magnetoacoustic waves using the 
parameters of MC. The analytical e xpressi ons of jwfll and \uji \ of some branches are present specifically. 

(a) Fast waves We see fr om F ig. 16(b)[ \Fx \ ^ \M\ (solid line) at large scales makes two fluids strongly coupled and oscillate 
together. The forces in Eq. dPli satisfy 

fi,x ^ fn^x ^ fi^z ^ fn,z (D3) 

Pi Pn Pi Pn 


to make neutrals move together with ions. The perpendicular component of force is larger due to the involvement of magnetic 
pressure. jF^I/l-^l (dashed line) increases towards/cdec.ni and exceeds |Pi|/|M| (open circles) at fcdec.ni- Within the cutoff interval, 
|Fa^|/|M| also increases to be larger than 1. The strong friction force acting on ions efficiently dissipates the oscillatory motions 
of ions. In the mean time, due to the weak coupling with ions, |P„|/|M| (filled circles) in neutrals substantially drops to the same 
level as |/|M| after the cutoff interval, and further decreases together with |P 2 |/|M| towards smaller scales. Accordingly, no 
wave motions can be developed in neutrals after they decouple from ions at fcdec.ni- 
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At scales smaller than 1 /fcdec, in, the forces have relations 


fi,x ^ fn,x fi,z ^ fn,z 

Pi Pn Pi Pn 

fi^x fi,z-! fn,x » fn,z- 


(D4a) 

(D4b) 


With decreasing frictions in both directions, ions oscillate under the pressure forces, and margina lly affe cted by neutrals. The 
damped fast waves reemerge, but only in ions, with a phase speed VAi instead of Va (shown in Fig. |16(a)| l. 

(b) Usual slow waves O wning t o the constant propagation direction {6 = 45°) adopted, we observe different wave behavior 
of the slow waves here (Fig. |16(c)| l from that in Fig. |7] despite the same environment parameters used. First, since the relation 
k\\ k± ^ k breaks down, fcdec.ni is smaller than (see Eq. (l52l l and (l5Tt ). Thus fcdec.ni becomes the cutoff scale of the 
usual slow waves in this case. Second, according to the expression in Table @1 the ratio fcc,ti/fcc,t 2 substantially increases due 
to the relatively large value of cos0. Meanwhile, the ratio |a;/|/|a;_R| of the ion slow waves within [fcc.ti, fcc,t 2 ] also increases, 
which leads to more efficient damping and a shorter wavenumber interval of this wave branch. In addition, with a constant cos 9, 
we cannot observe a decreasing \uji\ with k. From above comparison we clearly see that the wave propagation direction can 
significantly influence the slow wave behavior. 

Similar analysis to fast modes applies to the usual slow modes on scales larger than l/fcdec,ni- Within the interval [/cdec.ni, fcc.ti], 
ions are still strongly coupled with neutrals and undergo the same forces as neutrals, which we will discuss later for the neutral 
slow waves. The propagating ion slow waves arise at fcc.ti after |Pi|/|M| substantially i ncrease s and \M\ dominates over \Fx\. 
But due to the rise of iF’^l/IMj, the wave is strongly damped with |w/| > |w/j| (see Fig. |16(c)| l and cut off shortly at fccd 2 - The 
ion slow modes can only be fully resumed beyond fcdec.in when | Fz | /1M | decreases to 1. Then friction forces further decrease and 
become negligible compared with \M\ towards smaller scales. The forces have relations as 


fi,x ^ fi,z fuA ^ fn,x 

Pi Pi Pn Pn 


(D5) 


(c) Neutral slow waves As show n in F ig. |16(e) with 9 = 45 ° used, we can clearly see the change in the damping rate jw/l 
of the neutral slow modes (see Eq. ( I58bl i and (I59ll). Eig. |16(f)| illustrates the forces corresponding to the neutral slow modes. 
|P„|/|MI significantly increases and results in a nearly isotropic force on neutrals. The compression produced by the usual slow 
modes at larger scales drives oscillations in neutrals. Without much interference from ions, isotropically propagating sound waves 
are able to arise and maintain in the neutral fluid. Before reaching fcdec.in, ions are still coupled with neutrals. As a results, the 
forces are 


fn,x 

Pn Pn Pi Pi 


(D6) 


At scales smaller than k^^ | Fj | /1M | drops and converges with | Fz | /1M |. Accordingly fi^x!Pi further decreases and becomes 

significantly smaller compared with the forces on neutrals. 


THE EOKKER-PLANCK DIEEUSION COEEEICIENTS 

(a) Super-Alfvenic turbulence Since the hydrodynamic turbulence over [F, Ia\ does not contribute to particle scattering, the 
integral interval of is from the injection scale of GS95 turbulence Ia to the damping scale. In the following expressions of 
Dfifi for the case of super-Alfvenic turbulence, we use I a as the injection scale L and Ma = 1 applies correspondingly. 

In Alfvenic turbulence, by applying NET, we can get 




4z;vTr(l - p?) 
3FF2 



-ro/3 


dx 



/tVa 2 dr] Jl{w) 
-) — r - 


^11 rjxR^ 


(El) 


where ’G’ refers to gyroresonance. Here the dimensionless quantities are: R = u/(HF) (CR rigidity), x = kL, p = cos0. 
And = Av\\/v represents the dispersion in particle pitch angles due to the nonlinear effect. Jn{w) is Bessel function with 
w = kxv_L /il. In QLT, with the J-function resonance condition used, becomes 




X /’^1 , pVA^2d9 Jl{w) 

dx I (IH-)-^exp 
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(E2) 


Eor fast modes, we obtain 




FF2 




Ap^ 


(E3) 


^^;.(NLT) 
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U ni 


(a) Wave frequencies of fast modes 



kVA 


U ni 

(b) Forces of fast modes 



kVA 


y ni 



kVA 
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(c) Wave frequencies of usual slo\v modes 


(d) Forces of usual slow modes 



y ni 



y ni 


(e) Wave frequencies of neutral slow modes 


(f) Forces of neutral slow modes 


Fig. 16.— Wave frequencies and force ratios vs. normalized k of magnetoacoustic waves using MC pai'ameters. In (a), (c), and (e), the propagating |) and 
damping (| cjj |) components of wave frequencies are represented by solid and dashed lines respectively. In (b), (d), and (f), the displayed force ratios are | Fa; | /1M | 
(solid line), |F 2 |/|M| (dashed line), |Fj|/|M| (open circles), and |Fn,|/|Af| (filled circles). The horizontal dashed lines indicate the position |/|/|M| = 1. 


corresponding to gyroresonance and TTD interaction. For low-rigidity CRs with their tl shorter than the damping scale of 
turbulence, Jn{'^) can take the asymptotic form {wl2)'^ln\ at small argument. Then _Dj^(NLT) approximately becomes 
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The in QLT takes the form 


D^^iQLT) = 


Z)7r(l - 


LR^ 


X-^'^dx [\l + (x- 

Jl Jo VT] fj. \ 


rjfiR 


Its simplified form for low-rigidity CRs is 


D^:^(QLT) = 


VTry/JI{l — fi^) f2 2y^ 1 - /i2 
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ALy/R 


21/r2 


Diffusion coefficients in the case of slow modes are 
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(E7) 
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(b) Sub-Alfvenic turbulence Sub-Alfvenic turbulence have different energy spectra for weak and strong turbulence regimes. 
Accordingly, the diffusion coefficient is a sum of the components in both weak and strong turbulence. 

Eor the gyroresonance with Alfven modes, we have ^ g. Here 


^i.^(NLT) = 


2v^{l _2 tiVA.2 1 JlM 

■ / x_L dx±{l + - --^ 

Ji V A/t 


LR^ 


exp 
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(Ell) 


is D'J,, in weak turbulence, and 

/x/x 


ZJ« ,,(NLT) = 
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is in strong turbulence. We can see ZJ®^ j,(NLT) is similar to the ZJ®^(NLT) in super-Alfvenic case (Eq. (lEUl L but with a 
factor added (Ma < 1). 

In QLT, since weak turbulence only has smaller structures developed in perpendicular direction, i.e. fey = L~^ over [L, Itr], 
weak turbulence does not contribute to gyroresonance scattering unless CRs have comparable to L. Thus we only consider 
the contribution of strong turbulence and get 


^^^,.(QLT) = 


4x71(1 f\l I dVA^^dyJHw) 
iLR? Jo ^ V r]]i w'^ 


exp - 


X|| 
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rjfiR 


(E13) 

The energy spectrum of fast modes are not dependent on turbulence regime s. Ea st modes have similar diffusion coefficients in 
sub-Alfvenic turbulence as those in super-Alfvenic case. The expressions Eq. (lE31 l- (l74l i still apply, but here L is the driving scale 
of turbulence and should be multiplied accounting for the ratio between kinetic and magnetic energies of turbulence. 


Slow modes have 


ZJG ,„(NLT) = 
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in weak turbulence, and 
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in strong turbulence. 

QLT in strong turbulence gives 


i9;i(QLT) = 
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